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Annually, farmers in southern Africa manage their land resources and prepare their fields 
for cultivation by burning crop residual debris, with a peak in the burning season 
occurring during August and September. The emissions from these fires in southern 
Africa are among the greatest from fires worldwide, and the gases and aerosol particles 
produced adversely affect air quality large distances from their source regions, and can 
even be tracked in satellite imagery as they cross the Atlantic and Pacific Ocean basins. 
During August and September 2000 an international group of researchers participating in 
the Southern African Regional Science Initiate field experiment (SAFARI 2000) made 
extensive ground-based, airborne, and satellite measurements of these gases and aerosols 
in order to quantify their amounts and effects on Earth’s atmosphere. 

In this study we interpreted the measurements of smoke aerosol particles made during 
SAFARI 2000 in order to better represent these particles in a numerical model simulating 
their transport and fate. Typically, smoke aerosols emitted from fires are concentrated by 
mass in particles about 0.3 micrometers in diameter (1,000,000 micrometers = 1 meter, 
about 3 feet); for comparison, the thickness of a human hair is about 50 micrometers, 
almost 200 times as great. Because of the size of these particles, at the surface they can 
be easily inhaled into the lungs, and in high concentrations have deleterious health effects 
on humans. Additionally, these particles reflect and absorb sunlight, impacting both 
visibility and the balance of sunlight reaching Earth’s surface, and ultimately play a role 
in modulating Earth’s climate. 

Because of these important effects, it is important that numerical models used to estimate 
Earth’s climate response to changes in atmospheric composition accurately represent the 
quantity and evolution of smoke particles. In our model, called the Community Aerosol 
and Radiation Model for Atmospheres (CARMA) we used estimates of smoke emissions 
based on field studies and observations made with the NASA Terra and TRMM satellites. 
The meteorology used to calculate the transport was based on an assimilation of observed 
meteorological conditions provided by the National Center for Atmospheric Research. 
Our model was able to simulate much of the observed day-to-day variability in the smoke 
aerosol loading observed over the continent by the NASA AERONET network of 
ground-based sun photometers, suggesting that the observed variability is due more to 
meteorology than to day-to-day variability in emissions. We compared the simulated 
vertical profile of smoke concentrations with measurements made with NASA ground- 
based (MPLNet) and airborne lidars (Cloud Physics Lidar) and the NASA Ames 



Airborne Tracking Sun Photometer and found that to a good approximation the smoke 
aerosols are emitted in a well-mixed layer near the Earth’s surface called the planetary 
boundary layer. Because of the relative geographic sparseness of the AERONET, 
MPLNet, and airborne observations, it is important to also look at the observations made 
from satellites. Here we looked at observations of the smoke plumes made with the 
NASA MODIS, MISR, and TOMS satellite instruments, all of which provide important 
information about the quantity and distribution of aerosols present. These satellite 
instruments see the aerosols in somewhat complimentary ways, and we found important 
differences in their observed aerosol amounts, particularly over land, which highlights the 
difficulty in making these measurements. Overall, CARMA compared well with the 
observations available and we were able to constrain some of the parameters needed to 
accurately simulate the evolution of the optical properties of smoke aerosols as they are 
transported over long distances. 
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14 Abstract 

15 This study investigates the transport and optical properties of southern African biomass burning 

16 aerosols using an offline three-dimensional aerosol microphysical and transport model. Here we 

17 use Sun-photometer retrieved particle size distributions and monthly mean, satellite derived 

18 smoke emissions as input parameters in our model. We find that using these observations in our 

19 model allow us to reproduce the measured optical properties of smoke aerosols collected during 

20 the Southern African Regional Science Initiative campaign (SAFARI 2000). In particular, we 

21 find similar day-to-day oscillations in the simulated aerosol optical thickness (AOT) relative to 

22 AERONET measurements, suggesting that variations in aerosol loading are controlled more by 

23 transport processes than fluctuations in aerosol emissions. On a monthly basis, the model and 

24 observations from the MODIS, MISR, and EP-TOMS show that the dominant transport of smoke 

25 during September 2000 was westward towards the Atlantic Ocean. However, the observations 

26 show higher AOT values over the ocean than the model. These higher values observed by the 

27 satellites may be from the contribution of other aerosols or the condensation of gases onto the 

28 smoke, which are not simulated in this study. These higher values may also be a result of poor 

29 single scattering assumptions and/or the contamination of sub-pixel clouds in the satellite 

30 retrievals. At locations over Africa, we find discrepancies in AOT comparisons among the model 

31 and satellites which limits our ability in determining the performance of our model. We also 

32 suggest strategies for improving the treatment of particle size distributions of smoke aerosols in 

33 microphysical models. 

34 

35 
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36 Index Terms: 0305 Atmospheric Composition and Structure: Aerosols and Particles; 0312 

37 Atmospheric Composition and Structure: Air/sea constituent fluxes; 0322 Atmospheric 

38 Composition and Structure: Constituent sources and sinks; 0368 Atmospheric Composition and 

39 Structure: Troposphere - Constituent transport and chemistry; 9305 Geographic Location: 

40 Africa; Keywords: biomass burning smoke aerosols, transport model, biomass burning smoke 

41 aerosol sources, smoke removal mechanisms. 



42 1. Introduction 

43 Biomass burning is a significant regional and global source of aerosols and trace gases 

44 [ Crutzen and Andreae, 1990; Hao et al., 1990; Andreae, 1993]. The aerosols from both natural 

45 and anthropogenic burning of vegetation impact atmospheric chemistry, the radiation budget, air 

46 quality, and carbon cycling in ecosystems. It is now well known that biomass burning fires emit 

47 large amounts of carbon dioxide (CO2), carbon monoxide (CO), nitrogen oxides (NO x ), methane 

48 (CH 4 ), nonmethane hydrocarbons (NMHC), oxygenated volatile organic compounds (OVOC), 

49 and reduced nitrogen-containing species (i.e., ammonia (NH3) and hydrogen cyanide (HCN)) 

50 [Crutzen and Andreae , 1990 and Yokelson et al., 1996]. Some of these gases contribute to the 

5 1 greenhouse effect or are involved in photochemical processes that will generate greenhouse 

52 gases. In particular, biomass burning produces large amounts of CO (300 to 700 Tg(CO) y 1 ) and 

53 NO x (~8 Tg(N) yr~ l ), which are important precursors to the production of tropospheric ozone 

54 (O3) [ Seinfeld and Pandis, 1998]. Aerosols produced from biomass burning also influence the 

55 radiative properties of the atmosphere. Directly, biomass burning aerosols affect the radiation 

56 budget by scattering and absorbing shortwave (solar) radiation. Indirectly, these aerosols can 

57 serve as effective cloud condensation nuclei (CCN) at supersaturations greater than 0.5% which 

58 leads to the modification of cloud properties [Warner and Twomey, 1969]. IPCC [2001] 

59 provided a summary estimate that the global mean radiative forcing due to biomass burning 

60 aerosols is -0.2 W nf 2 with an uncertainty as large as a factor of 3 and a “very low” level of 

61 scientific understanding. Additionally, satellites often observe smoke plumes being transported 

62 over intercontinental and continental scales, thereby degrading visibility not only at the smoke 

63 source but also downwind of the smoke source [Posfai et al., 2003; Stein et al., 2003]. In 

64 particular, the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) observed a significant smoke 
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65 event, no\v referred to as the “River of Smoke,” [Annegarn et al. 2002] where aerosols from 

66 biomass burning was transported over southern Africa and the Indian Ocean on 4 September 

67 2000 as illustrated in Figure la. 

68 Central and southern Africa experiences the most extensive biomass burning in the world 

69 [ Scholes & Andreae, 2000]. It is estimated that approximately 49% of atmospheric carbon 

70 produced by fires worldwide comes from the savanna fires of Africa [Dwyer et al, 2000; Scholes 

71 & Andreae, 2000] with emissions ranging from 80 to 785 Tg(C) yr _1 [ Scholes et ah, 1996, 

72 Barbosa et al., 1999; van der Werfet al., 2003]. There have been substantial scientific efforts at 

73 estimating smoke emissions, however the components (i.e., burned area, fuel loading, and 

74 combustion factors) used to calculate the emissions are highly uncertain producing values that 

75 fall within a large range. 

76 The Southern African Regional Science Initiative campaign (SAFARI 2000) was an 

77 international scientific campaign to study land-atmosphere interactions in southern Africa during 

78 the wet season in 1999 and the dry season in 2000 [Swap et al., 2003]. One of the main 

79 objectives of SAFARI 2000 was the characterization and quantification of regional emission 

80 sources from savanna burning using surface, aircraft, and remote sensing measurements of 

81 aerosol chemical composition, vertical distributions, refractive indices and particle size 

82 distributions. During the early part of the dry season in August, these fires were prescribed as a 

83 land management tool. The burning of wet soils at the beginning of the dry season produces low 

84 fire intensities, which result in limited vegetation consumption and limited damage to the soil. 

85 Farmers bum extensively in the early dry season to achieve rapid nutrient release prior to the 

86 new growing season and to stimulate regrowth of palatable grasses for their cattle. Early burning 

87 is also used in national parks as a preventative measure against late dry season fires, which tend 
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88 to be more intense and destructive. The estimated mean fire size, maximum fire size, and total 

89 burned area over Africa during the SAFARI 2000 campaign were 7.83 km , 1458 km , and 

90 95,962 km 2 , respectively [Hely et al., 2003]. 

91 This paper focuses on modeling the transport and optical properties of biomass burning 

92 aerosols over continental Africa during SAFARI 2000. The optical properties investigated here 

93 include aerosol optical thickness (AOT), Angstrom exponent (a), single scattering albedo (SSA) 

94 and vertical extinction profiles. Our first goal in this study is to determine suitable assumptions 

95 about aerosol emissions, injection altitudes, and initial particle size distributions that reproduce 

96 the observed optical properties of biomass burning aerosols. Our second goal is to better 

97 understand the microphysics needed to accurately simulate biomass burning aerosol properties. 

98 We conducted several simulations to investigate the model’s sensitivity to: the timing and 

99 altitude of the smoke emissions; the transport of a particle size distribution; microphysical 

100 processes; removal schemes; and smoke optics. These sensitivity tests will aid us in 

101 understanding the parameters that control our model results. 

102 Below, section 2 describes our model and input parameters used to simulate the biomass 

103 burning aerosols in our base simulation. Section 3 presents the results from our base simulation 

104 and sensitivity tests, and compares the results to observations made during SAFARI 2000. Lastly 

105 in section 4, we conclude with a discussion of the current ability to successfully model biomass 

106 burning aerosols produced from savanna fires over continental Africa. 

107 

108 2. Model Description 

109 2.1 Dynamical Component 
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110 Opr model is driven by meteorological fields derived from the National Center for 

111 Atmospheric Research (NCAR) Model for Atmospheric Transport and Chemistry (MATCH) 

112 [Rasch et al., 1997]. MATCH is an offline chemical transport model that uses wind, 

113 temperature, and pressure fields provided by the National Center for Environmental Prediction 

1 14 (NCEP) reanalysis package (NCEP/NCAR reanalyses [Kalnay et al, 1996]) and the physical 

115 parameterizations from the NCAR Community Climate Model version 3 (CCM3) [Kiehl et al., 

116 1996]. The NCEP/NCAR reanalyses are gridded at T63 horizontal resolution (approximately 

117 1.875° x 1.875°) with 28 vertical sigma layers extending from the surface to approximately 35 

118 km, and are available each day at 0000, 0600, 1200, and 1800 UTC. MATCH is run with half- 

119 hour timesteps, and the input fields are linearly interpolated to the current timestep. At each 

120 timestep MATCH employs the CCM3 physics to diagnose fields required to compute planetary 

121 boundary layer transport, convective mixing, and cloud and precipitation fields. We archive 

122 these fields every six hours for use in our aerosol microphysics and transport model. More details 

123 of our use of the MATCH derived fields are discussed in Colarco et al. [2003a]. 

124 

125 2.2 Aerosol Transport Component 

126 The global dynamical fields generated by MATCH are imported into our aerosol 

127 microphysical, transport, and radiation model, which is a version of the Community Aerosol and 

128 Radiation Model for Atmospheres (CARMA) developed at NASA Ames Research Center and 

129 the University of Colorado [Toon et al., 1988; Ackerman et al., 1995]. CARMA and MATCH 

130 have been previously applied to three-dimensional simulations of the transport and evolution of 

131 dust and carbonaceous aerosols [Westphal et al. [1991]; Colarco et al., 2002; 2003a, b; 2004]. 
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132 CARMA is a bin-resolving microphysical model that solves the aerosol continuity 

133 equation for source, transport, removal and transformation processes. The particle size 

134 distribution in CARMA is treated using a number of discrete bins distributed in radius space. In 

135 this study, the injected aerosols are distributed across 16 size bins spaced logarithmically 

136 between 0.01 jim and 10 pm radius sizes. For the smoke injection, we use a bimodal lognormal 

137 distribution based on observations (see section 2.6 and Figure 2). Additionally, only smoke 

138 aerosols are transported in the model (i.e., no dust or sea salt aerosols). The model can include all 

139 basic microphysical processes affecting aerosols; however in this study we only consider 

140 advection by winds, sedimentation, dry deposition, wet removal, coagulation, and prescribed 

141 aerosol sources. These processes are described in Colarco et al. [2003a]. Briefly, the advective 

142 and diffusive transports are solved using a mixed form of a numerically accurate, highly non- 

143 diffusive piecewise polynomial algorithm [Collela and Woodward , 1984] and a semi-lagrangian 

144 transport step following Lin and Rood [1996]. Particle transport by sedimentation and dry 

145 deposition are incorporated into the vertical transport component through fall velocities and 

146 deposition velocities, respectively. The fall velocity is computed for each size bin following the 

147 treatment by Pruppacher and Klett [1997]. The surface layer dry deposition velocity is 

148 calculated with a two-layer dry deposition model which accounts for sedimentation, molecular 

149 diffusion, and turbulent deposition across the lowest model layer [ Shao , 2000]. CARMA’ s size- 

150 dependent, wet removal (scavenging by precipitation) parameterization uses the cloud and 

151 precipitation fields calculated in MATCH and is treated as a first order loss process as described 

152 by Barth et al. [2000]. The coagulation of the smoke aerosols is treated via Brownian motion. 

153 Here, our coagulation algorithm preserves the total aerosol volume but decreases the number 

154 concentration of particles as small particles stick together. In general, coagulation modifies the 
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155 particle size distribution by increasing the mean radius of the particle size distribution at a rate 

156 approximately proportional to the air temperature, the square of the particle number 

157 concentration, and the inverse of the particle radius. 

158 For this study, CARMA is run at the same temporal and spatial resolution as MATCH, 

159 but on a limited area grid (180°W to 180°E and approximately 15°N to 60°S). The model is run 

160 for a number of simulations in order to test sensitivity to various parameterizations of sources 

161 and transport processes. The input smoke sources are calculated off-line, and are described in the 

162 next section. CARMA is driven for the period of August and September 2000 with half-hour 

163 timesteps. Here, the 6-hr input fields from MATCH are linearly interpolated to the desired 

164 timestep in CARMA. We include August to allow a 30-day spin-up time for aerosols in the 

165 model, and analyze only the September results. 

166 

167 2.3 Smoke Source Component 

168 In this study, we are interested in testing how our modeled smoke fields compare to 

169 observations relatively near fire sources. The aerosol emissions used in the model were 

170 interpolated from the Global Fire Emissions Database (GFED), which provides a global 1° x 1° 

171 gridded monthly mean map of fire carbon consumption for the period January 1997 to December 

172 2002 [van der Werfet al. , 2003]. The carbon consumed was constructed using the Tropical 

173 Rainfall Measuring Mission (TRMM) satellite and Visible Infrared Scanner (VIRS) hot spot data 

174 [Giglio et al., 2002], Moderate Resolution Imaging Spectroradiometer (MODIS) burned area 

175 information [Kaufman et al., 2003], and a spatially dependent relationship between fire counts 

176 and biomass consumed developed from the Camegie-Ames-Stanford Approach (CASA) 

177 biogeochemical model. To convert carbon consumed due to fires into aerosol emissions for our 
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178 model we assumed an emission factor of 10.0 g Total Particulate Matter (TPM) per kg of Dry , 

179 Mass (DM) [Sinha et al., 2003], which is representative of savanna and grassland vegetation in 

180 Africa assuming 45% of the DM burned was carbon [ Andreae and Merlet, 2001], The smoke 

181 emissions are sensitive to the emission factor, and there are uncertainties in the GFED dataset 

182 particularly with relation to the burned area estimates [van der Werfet al ., 2003]. However, 

183 sensitivity simulations with other emission datasets (e.g., Korontzi et al. [2003], not shown) yield 

184 qualitatively similar spatial distributions of the smoke sources. We also assumed the fires emitted 

185 aerosol at the same rate each day. Over the approximately hemispheric domain of the model, the 

186 total aerosol emissions are about 4.52 Tg in August and 7.92 Tg in September, where 

187 approximately 2.14 Tg and 3.51 Tg of the total emissions were produced over Africa in August 

188 and September, respectively. Figure 3 shows the spatial distribution of aerosol emissions 

189 interpolated to our model grid over Africa for September 2000, and the African locations 

190 analyzed in this paper. 

191 

192 2.4 Smoke Optical Property Calculations 

193 We calculated the AOT, SSA, Angstrom exponent, and vertical extinction profile of our 

194 simulated smoke aerosol distributions using Mie theory [ Wiscombe , 1979]. We represent the 

195 aerosols as spherical “smoke” particles, making no assumptions about their composition except 

196 their density (1,35 g cm"" 3 [Reid and Hobbs, 1998]) and their refractive index. We are not aware 

197 of any data suggesting that the SAFARI aerosols were external mixtures, so our approach is 

198 consistent with treating the aerosols as internal mixtures of absorbing and non-absorbing 

199 components. For simplicity, we selected a wavelength independent refractive index. The 

200 refractive index selected for the optical property calculations represents the mean complex 
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201 refractive index of days dominated by smoke aerosols at Ndola, Zambia for the month of 

202 September 2000 as retrieved by the Aerosol Robotic Network (AERONET) [Holben et al., 1998] 

203 Sun-photometers (see Table 1 and Section 3.1.1 below). Our criterion for classifying smoke 

204 dominated days was based on an analysis of the AERONET climatology of cloud-screened, 

205 quality assured aerosol optical thickness for southern African observations. Here we found that 

206 the monthly mean AOT at 500 nm [AOT 500 ] exceeded 0.3 during months dominated by biomass 

207 burning emissions. Accordingly, when the daily average of AERONET observations show 

208 AOT500 greater than 0.3, the day is classified as smoke domi n ated. Additionally, the 

209 humidification of the smoke aerosols was treated in the optical property calculations following a 

210 parameterization described by Magi and Hobbs [2003]. 

211 

212 2.5 Base Model Description 

213 The input parameters used in our base model are summarized in Table 1 . In our base 

214 simulation, we assume the smoke particles have a lognormal size distribution using number 

215 fractions (Nf), geometric median radii On), and standard deviations ( 0 ) retrieved by an 

216 AERONET almucantar scan [Dubovik et al, 2000] at Ndola, Zambia on 16 September 2000 (see 

217 location in Figure 3 and details in Table 1). We selected this particular particle size distribution 

218 because of the site’s proximity to emission sources and because it has the highest daily mean 

219 AOT 500 of all the study sites in the smoke region during September 2000. This high optical depth 

220 suggests that it was dominated by locally generated smoke aerosols. Figure 2 (solid black line) 

221 presents the vertically column-integrated volume size distribution [pm 3 /ji.m 2 ] used in the base 

222 model. It is not clear that the AERONET particle size distribution selected here represents young 

223 smoke (i.e., less than 3 minutes old [ Haywood et al., 2003]). Rather, it is likely more 
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224 representative of a mixture of young and moderately aged smoke, as well as contributions from. 

225 other aerosols (e.g., dust), existing in different layers in the atmosphere. We find that a volume 

226 size distribution fit to parameters representative of young smoke (. Haywood et al. [2003], grey 

227 dashed line in Figure 2) has a narrower fine mode volume size distribution than the AERONET 

228 retrieval, which suggests that the size distribution selected for our base model may not be 

229 representative of only young smoke. On the other hand, the coarseness of our model’s spatial and 

230 temporal resolution prohibits a detailed treatment of truly fresh smoke aerosols. Therefore, the 

23 1 proximity of the Ndola site to emission regions suggests that this site provides a reasonable 

232 compromise between young and aged hazes at a scale our model can reasonably simulate. 

233 Eck et al. [2003] observed significant diurnal variability in AOT 500 over six sites in 

234 Zambia (Mongu, Ndola, Senanga, Mwinilunga, Zambezi, Solwezi) with a minimum AOT 500 at 

235 0900 UTC and maximum at 1500 UTC. For southern Africa, the TRMM satellite, which makes 

236 two overpasses every two days at the equator, also indicated a strong diurnal cycle in fire counts 

237 in August and September with a maximum number of fires occurring approximately 3 hours after 

238 local solar noon [ Giglio et al., 2003]. To represent a diurnal cycle in our base model, the aerosol 

239 emissions are concentrated during a 12-hour period between 0600 and 1800 UTC. Additionally, 

240 the smoke emissions are mixed throughout the depth of the planetary boundary layer to account 

241 for the convection and buoyancy associated with the heat and moisture in the fires [ Liousse et al., 

242 1996; Lavoue et al., 2000]. The smoke emissions were evenly distributed between the surface 

243 model layer and model layer associated with the top of the planetary boundary layer (PBL) 

244 reported by MATCH. On average, the MATCH planetary boundary layer height (PBLH) at the 

245 six study sites was lowest between 0000 and 0600 UTC, with heights at approximately 0.3 km, 
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246 and peaked at 1200 IJTC. The peak height of the MATCH PBL ranged between 1 .5 km and 3.7 

247 km. 

248 

249 3. Model Results and Discussion - , 

250 The simulated optical properties are compared to satellite, aircraft and ground-based 

25 1 measurements made over southern Africa during September 2000. For this study, the model 

252 results are compared to data at six locations: Etosha Pan, Namibia; Inhaca, Mozambique; 

253 Mongu, Zambia; Ndola, Zambia; Senanga, Zambia; and Skukuza, South Africa (Figure 3). 

254 These locations were selected based on the amount of available data obtained during the 

255 SAFARI 2000 campaign and their proximity to aerosol source regions. 

256 In section 3.1, we compare the results from the base model to observations made during 

257 SAFARI 2000, and discuss how well the initial input parameters used in our base model 

258 represent the aerosol distributions near the smoke source regions. Section 3.2 presents the results 

259 from a series of sensitivity tests to understand the influence the smoke source function and 

260 microphysical processes have on the model results. 

261 

262 3.1 Results - Base Model 

263 3.1.1 Comparisons to Ground-Based Optical Measurements from AERONET 

264 During SAFARI 2000, there were 17 Sun-sky photometer sites active in southern Africa 

265 as part of the globally distributed, federated Aerosol Robotic Network (AERONET) [Holben et 

266 al., 1998], The standardized AERONET instrumentation measures the direct solar beam 

267 transmission to determine the column-integrated AOT, and additionally retrieves water vapor, 

268 aerosol particle size distributions, and aerosol single-scattering albedo [ Holben et al, 1998; 
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269 Dubovik et al., 2000]. In this study, we are using level 2.0 cloud-screened quality-assured data 

270 [Smirnov et al, 2000]. The estimated uncertainty in AOT ranges between +0.010 and +0.021, 

271 which is spectrally dependent with higher errors in the UV [ Eck et al., 1999]. The uncertainty in 

272 SSA is estimated to be +0.03. The expected accuracy in the column-integrated volume size 

273 distributions are approximately 15-25% for radius sizes between 0.1 pm and 7 pm and 25-100% 

274 (or less than 10% of dV/d( In r) in maximum) for sizes less than 0.1 pm and greater than 7 pm 

275 [Eck et al., 1999; Dubovik et al., 2000], 

276 Figure 4 presents AERONET and simulated AOT 50 o for each day in September 2000 at 

277 the six African study sites noted above. We point out that the AERONET measurements are 

278 essentially point measurements, while the resolution of our model grid box is approximately 200 

279 km x 200 km. To make the model and measurements more comparable, we show all the 

280 retrievals made by AERONET (asterisks) in a single day and the daily mean AOT 500 simulated 

28 1 by the model. The solid line represents the results from the base model, and the dashed line 

282 represents the results from a simulation where we generated variable emissions on a daily basis 

283 (see Section 3.2). Qualitatively, we find that simulated AOT 500 falls within the daily range of 

284 AERONET retrievals on most days, and illustrates a similar temporal variation. The day-to-day 

285 fluctuations in the simulated AOT 50 o are remarkable given the fact that the emissions used in the 

286 model are released at a constant rate for the month of September. This suggests that the 

287 oscillations in the simulated AOT 500 are controlled more by the variability in the meteorological 

288 conditions then the temporal and spatial variations in the fire activity. This has also been 

289 observed by Myhre et al. [2003]. On days when the model does not agree with AERONET, the 

290 model is generally underestimating the daily mean AOT 500 . However, at Etosha Pan we do find 

291 that the model overestimates the daily mean in the latter half of the month. 
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292 We also find that the timing and magnitude in the simulated AOT 500 compares better to 

293 AERONET at Etosha Pan, Inhaca, Ndola, and Skukuza than at Mongu and Senanga. Although 

294 the model AOT 50 o does not correlate well with the AERONET observations at Mongu and 

295 Senanga, it is interesting to note that the model AOT 500 at Mongu is well correlated with the 

296 model AOT at Senanga. This is not necessarily surprising since the two sites are within a couple 

297 of hundred kilometers of each other (i.e., adjacent grid cells). However, the fact that the 

298 AERONET observations at the two sites are also well correlated with each other suggests that 

299 the aerosol loading at both locations is affected by dynamical and emission processes happening 

300 on synoptic rather than local scales. MOD IS (AOT 550 ) retrievals from 4-6 September and 15 - 

301 17 September (not shown) confirm the large-scale aerosol loading, where it shows AOT550 

302 values greater than 1.0, presumably associated with biomass burning, over a large portion of 

303 Zambia which encompasses both of these sites. Comparing the MODIS observations on these 

304 days to the model results (not shown) reveals that the model has the plume in the correct location 

305 during the 4-6 September but underestimates the AOT relative to MODIS. However during the 

306 15 - 17 September, the model plume is placed too far south relative to MODIS, thereby causing 

307 poor correlations between the model and AERONET during this period. These results suggest a 

308 problem with the model dynamics during this period, although it does not rule out that there may 

309 be a major and perhaps episodic aerosol source not in the GFED dataset. 

310 Figures 5-7 present scatterplots of simulated and AERONET daily mean AOT 500 , 

311 SSA440, and Angstrom exponent [CX440/870] on smoke dominated days at all six study sites, 

312 respectively. The correlation coefficient (R) and the equation of the line of best fit are also 

313 reported in each scatterplot. Although the model AOT is well correlated with the AERONET 

314 observations (R=0.78, Figure 5), we find that the model underestimates the daily AOT 500 by up 
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315 to a factor of 2 when AERONET retrievals are greater than 1.0. This may be a result of spatially 

316 small, dense aerosol plumes observed by AERONET which cannot be resolved in our model. 

3 17 When we discard the AERONET observations greater than 1 .0, we find that the correlation 

318 coefficient increases to 0.81 and the slope of line of best fit is closer to 1 (slope = 0.80). In the 

319 comparisons of daily SSA 440 and Angstrom exponent, we find that the model results do not show 

320 as much variability as the AERONET retrievals (Figures 6 and 7, respectively). Here the 

321 simulated SSA 440 ranges between 0.85 and 0.90, where, the AERONET retrievals range between 

322 0.85 and 0.95. Given the large uncertainty in the AERONET SSA retrievals (i.e., +0.03 and 

323 represented by single bar at top-left comer of Figure 6), the model is within the measurement 

324 error for SSA. The simulated Angstrom exponent ranges between 1 .70 and 1.95, where the 

325 AERONET retrievals range between 1.40 and 2.10. This difference suggests that our model is 

326 not capturing the variability in the size distribution correctly. The model results may be 

327 improved if we consider the condensation of gases onto the smoke aerosols, and the contribution 

328 from other aerosol species that exist at these study sites (i.e., Aeolian dust, sea salt, and fossil 

329 fuel combustion). The model’s sensitivity to initial size distribution will be discussed further in 

330 section 3.2. 

331 Figures 8 and 9 compare the monthly mean spectral dependence of AOT and SSA (i.e., 

332 how AOT and SSA changes with wavelength) of smoke dominated days for September 2000, 

333 respectively. The asterisks connected with a dashed line represent the AERONET data and the 

334 solid line represents the results from the base model. For the remaining portion of this paper, we 

335 only discuss the results from Ndola and Inhaca since we find similar results at the other sites. 

336 Etosha Pan, Mongu, and Senanga comparisons are similar to the results found at Ndola, while 

337 the Skukuza comparison is similar to the results found at Inhaca. At Inhaca and Ndola, we find 


16 



338 that the mean wavelength dependence of AOT simulated by the base model is similar to the 

339 wavelength dependence retrieved by AERONET. The maximum difference between the model 

340 and AERONET monthly mean AOT at all measured wavelengths is approximately 30% at both 

341 sites. We find that the simulated wavelength dependence of SSA is comparable to AERONET at 

342 Ndola but not at Inhaca. At Inhaca, we find that the simulated monthly mean SSA is within the 

343 uncertainty of the AERONET retrievals at wavelengths less than 870 nm. However, at 

344 wavelengths greater than 870 nm, the model overestimates the absorption by approximately 

345 50%. Again, this suggests that our model is not capturing the variability in the size distribution 

346 and/or refractive index correctly. Also note that the SSA is wavelength dependent. Here we find 

347 that it is larger at shorter wavelengths for both the model and the data despite the refractive index 

348 for the model being assumed to be wavelength independent. Further investigations of wavelength 

349 independent and dependent refractive indices are discussed in section 3.2. 

350 The spectral dependence of the AOT lends important insight into the size distribution of 

351 particles being considered. Using AERONET data, Eck et al. [1999] characterized the spectral 

352 properties of African and South American biomass burning aerosols. A typical method of 

353 inferring particle size information from Sun-photometer observations is to compute the 

354 Angstrom exponent (a), which requires the AOT retrieval at two wavelengths (see Eck et al. 

355 [1999], equation 2) and represents the first derivative of the variability of log e (AOT) versus 

356 loge(wavelength). As defined, Angstrom exponent values near 2.0 correspond to measurements 

357 that are dominated by small-sized particles (i.e., smoke aerosols), while values near 1.0 

358 correspond to large-sized particles (i.e., dust and sea salt). Since the major contribution to AOT 

359 for biomass burning is from sub-micron, accumulation mode size particles, Eck et al. [1999] 


360 recommend that the Angstrom exponent be computed from shorter wavelength pairs which have 
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361 much greater sensitivity to accumulation mode size particles than from longer wavelength pairs-. 

362 To acquire information about the relative contributions of the accumulation and coarse mode 

363 particles to the total AOT, Eck et al. [1999] also reco m mend computing a curvature term (a’) 

364 which provides information on the second derivative of the log e (AOT) versus log e (wavelength) 

365 relationship. As defined in Eck et al. [1999], equation 8, larger (positive) values of the curvature 

366 term indicate relatively larger contributions of accumulation mode versus coarse mode particles 

367 to the total AOT. For biomass burning aerosols in Mongu, Zambia, Eck et al. [1999] find a = 

368 1.21 for the 380 nm and 440 nm wavelength pair and a curvature value of a’ = 1 .24. 

369 Table 2 lists the monthly mean AOT 50 o, oc 38 omo> ot’, SSA 440 , SSA 670 , SSA 870 , and SSA 1020 

370 of smoke dominated days at Inhaca and Ndola. To compute the derivative of the Angstrom 

371 exponent (a’), we used wavelengths at 380, 500, and 870 nm. The AERONET observations at 

372 both locations show similar Angstrom exponent values, which for this short-wavelength pair 

373 indicate similar fine mode sizes, but very different curvature values. The slightly lower 

374 Angstrom exponent and much lower curvature value at Inhaca relative to Ndola is consistent 

375 with the known contribution of coarse mode aerosols at this location [Eck et al., 1999], while the 

376 Ndola observations are more consistent with a biomass burning dominated site. The base model 

377 results at both locations are similar to each other, and the model agrees reasonably well with the 

378 observations at Ndola, but underestimates the coarse mode contribution to the AOT at Inhaca. 

379 The fact that the model underestimates the coarse mode contribution to AOT at Inhaca is not 

380 surprising since we are not simulating the sea salt and dust aerosols, which is likely to explain 

381 the observations. As shown in Figure 9 and reported in Table 2, the model is slightly less 

382 absorbing at Ndola and slightly more absorbing at Inhaca compared to AERONET. However, the 

383 model SSA is within the uncertainty of the AERONET retrieval error bars at all wavelength 
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384 retrievals, except for the 1020 nm wavelength retrieval at Inhaca. The 1020 nm retrieval may be 

385 understood as due to the coarse mode aerosols explained above. 

386 

387 3.1.2 Comparisons of Vertical Extinction Profiles from Airborne and Ground-Based 

388 Measurements ■' 

389 During SAFARI 2000, the NASA 14-channel Ames Airborne Tracking Sun-photometer 

390 (AATS-14) was operated aboard the University of Washington’s Convair-580 research aircraft 

391 [ Schmid et al., 2003]. The AATS-14 made measurements of smoke originating from flaming 

392 grass fires occurring at several locations in central and southern Africa. The AATS-14 measured 

393 the transmission of the direct solar beam for 14 discrete wavelengths from 354 nm to 1557 nm 

394 from which the spectral AOT can be derived. In general, the AATS-14 profiles showed that most 

395 or all of the aerosols were below 4 km with extinctions at 525 nm as large as 0.35 km -1 during 

396 the SAFRI 2000 campaign [Schmid et al., 2003]. Depending on synoptic conditions, some 

397 profiles showed multiple-layers, where an aerosol layer would be present between the surface 

398 and approximately 1.5 km and an additional layer between 2 km and 4 km [Schmid et al., 2003]. 

399 The Cloud Physics Lidar (CPL) also made measurements during the SAFARI 2000 

400 campaign. The CPL was built for use on the NASA ER-2 high-altitude aircraft. This instrument 

401 provides information on cloud height and structure as well as aerosol optical thickness, and 

402 operates simultaneously at 3 wavelengths: 355 nm, 532 nm, and 1064 nm [McGill et al., 2003]. 

403 The CPL fundamentally measures range-resolved profiles of volume 180-degree backscatter 

404 coefficients, with a vertical resolution fixed at 30 m. During SAFARI 2000, the CPL profiles 

405 showed that most of the aerosols were below 4 km with extinctions at 532 nm as large as 0.8 

406 km -1 [McGill et al, 2003] . 
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407 The Micro-Pulse Lidar (MPL) was also operational during the SAFARI 2000 field 

408 campaign [ Spinhirne et ah, 1995]. The MPL is a single channel (523 nm), autonomous, eye-safe 

409 lidar system originally developed at the NASA Goddard Space Flight Center. Several 

410 instruments are collocated with AERONET Sun-sky radiometers as part of a global lidar network 

411 (MPLNET, Welton et al., [2001]). The MPL is used to determine the vertical structure of clouds 

412 and aerosols to produce optical properties, such as extinction and optical depth profiles of the 

413 clouds and aerosols. The MPL profiles at multiple African locations generally showed most of 

414 the aerosols below 4 km with extinctions at 523 nm as large as 0.3 km -1 . 

415 Figure 10 compares the simulated mean vertical extinction profiles at 525 nm to 

416 measurements collected by the A ATS- 14 and CPL at comparable wavelengths (i.e., CPL at 532 

417 nm) at Mongu and Senanga, Zambia on 6 September 2000. In this figure, the simulated mean 

418 extinction profiles represent the average model results at 0600 and 1200 UTC, the CPL profiles 

419 represent the average of four samples taken between 0800 and 1 100 UTC (black lines), and the 

420 AATS-14 profiles represent a single sample (red lines). The blue lines are the results from the 

421 base model, and the green lines are results from a sensitivity study testing different injection 

422 height strategies (see Section 3.2). In general, we find that the results from the base model 

423 compare well to the observations at the selected African locations. Again, we find similarities in 

424 the simulated and observed vertical profiles at both sites which may be associated with a large- 

425 scale smoke event. However, the AATS-14 measurement observes a well defined smoke layer 

426 between 2 km and 4 km where the base model does not. This may be a result of excess mixing in 

427 our model, which is possibly related to the excessively coarse vertical resolution in the model. It 

428 is however difficult to draw conclusions about the behavior of the model from single smoke 
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429 events suph as these aircraft samples since they may have detected isolated plumes on a scale 

430 much smaller than the model resolution. 

431 Figure 1 1 compares the simulated mean vertical extinction profile at 525 nm to the mean 

432 MPL measurements at 523 nm at Skukuza, South Africa. We are using level 1.5a data products 

433 from MPL (real-time cloud and aerosol data). The computed mean includes all days when MPL 

434 was operational (i.e., 1-3, 6, 7, 9, 10, 13, 14, 17 September 2000). The dashed line represents the 

435 MPL measurements and the remaining lines represent the results from various model 

436 simulations. In this section, we are only concerned with the solid black line which represents the 

437 results from the base model. Here we find that the base model results compare very well to the 

438 MPL measurements. 

439 

440 3.1.3 Comparisons to Spaceborne Measurements of Aerosol Optical Thickness - 

441 Figure lb shows the simulated daily mean AOT at 550 nm on 4 September 2000, the 

442 same day as the SeaWiFS satellite image shown in Figure la. Here we find that the model 

443 produces a similar pronounced aerosol transport corridor across southeastern Africa and over the 

444 Indian Ocean or the “River of Smoke”. Qualitatively, this simple comparison demonstrates our 

445 model’s capability of transporting smoke aerosols properly on a single day. 

446 Figures 12 - 14 compare the simulated monthly mean AOT from the base model to the 

447 retrieved AOT from polar orbiting satellite platforms: the Moderate Resolution Imaging 

448 Spectroradiometer (MODIS) [Kaufman et ah, 1997; Tanre et al., 1997], the Multiangle Imaging 

449 Spectroradiometer (MISR) [Diner et al., 1998; Martonchik et al, 1998, 2002], and the Total 

450 Ozone Mapping Spectrometer onboard the Earth Probe satellite (EP-TOMS) [Torres et al., 

451 1998], respectively, over southern Africa for September 2000. MODIS and MISR are both flying 
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452 on the NASA Terra spacecraft, with an equator crossing time of approximately 10:30 AM and ■ 

453 10:30 PM local time. Both instruments therefore have similar temporal coverage, although their 

454 fields of view are somewhat different, reflecting the different capabilities of each instrument. 

455 The EP-TOMS instrument has equator crossing times of about 10:00 AM and 10:00 PM local 

456 time. For all three instruments only the day-lit observations are useful for obtaining AOT 

457 retrievals. From theoretical sensitivity studies, the estimated uncertainty in the AOT for MODIS 

458 is the larger of +0.05 and [+ O.2AOT550] over land and +0.03 and [+ O.O5AOT550] over ocean 

459 [ Kaufinan et ah, 1997; Tanre et al, 1997], for MISR it is the larger of +0.05 and [+0.2AOTx«?] 

460 over the land and ocean [ Kahn et al, 2001], and for EP-TOMS it is the larger of +0.1 and 

461 [+O.3AOT500] over land and ocean [Torres et al., 1998]. 

462 In order to quantitatively compare the satellite retrievals to our model we developed an 

463 ad hoc aggregation and sampling strategy for computing the monthly mean of the satellite 

464 retrieved AOT. The retrievals are typically available at a much higher spatial resolution than our 

465 model: 10 x 10 km 2 for MODIS, 17 x 17 km 2 for MISR, and approximately 24 x 24 km 2 for EP- 

466 TOMS. These are commonly referenced as “Level 2” products. The various satellite retrievals 

467 are also generally available in a gridded (“Level 3”) product which typically has a 1° x 1° spatial 

468 resolution, but the quantitative validity of these products when used to construct a monthly mean 

469 can be slightly dubious in situations where only a few number of pixels (others being cloud- or 

470 glint-obscured) dominate the monthly mean. For this reason, our goal in this section is to 

471 develop a method to quantitatively compare the satellite retrievals to our model and the satellites 

472 to each other by constructing our own gridded product from the Level 2 retrievals. 

473 Our approach was to first aggregate the retrievals at their native resolution to our model 

474 grid for each day, eliminating all but the highest quality pixels (based on quality assurance flags 
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475 or other information associated with the retrieval product). This way we find the average 

476 retrieved AOT in each model grid box on each day, as well as the number of retrievals available 

477 and the standard deviation about the mean. The monthly mean at each grid box was constructed 

478 from the daily averaged grid box values, but only using days which exceeded some threshold 

479 number of valid retrievals (the pixel threshold, eliminating grid boxes which are dominated by 

480 only a few retrievals). The pixel threshold value is. a function of the satellite instrument (i.e., we 

481 require more pixels for the high resolution MODIS and fewer for the low resolution TOMS) and 

482 the satellite repeat cycle (near daily coverage for MODIS versus the 7 - 8 day repeat cycle for 

483 MISR). The monthly mean at each grid box is then the mean of the daily values passing the pixel 

484 threshold test, weighted by the number of pixels per day. This strategy is applied in an ad hoc 

485 manner in the sense that we selected each satellite’s pixel threshold by iteration, selecting where 

486 small variations in the pixel threshold did not lead to significant changes in our computed 

487 monthly mean. Accordingly, the pixel threshold for each sensor is: 15 for MODIS, 5 for MISR, 

488 and 3 for EP-TOMS. Finally, the model AOT was sampled near the satellite overpass time and 

489 only where the satellite daily grid box averages passed the pixel threshold test. The model 

490 monthly mean was then constructed by weighting it in the same fashion as the satellite monthly 

491 mean. Thus, for each satellite we have a unique monthly mean of the model AOT which can be 

492 compared to the satellite retrieval at the model resolution. 

493 To assess our sampling and aggregation method, we compared our MODIS results to the 

494 MODIS Online Visualization and Analysis System (MOVAS) product (not shown) 

495 (http://g0dup05u.ecs.nasa.gov/Giovanni/) and found that the patterns in the monthly mean AOT 

496 are qualitatively consistent between the two products. While we have not fully evaluated this 

497 technique, it addresses the non-trivial problem of scaling and aggregating high resolution 
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498 satellite observations for quantitative comparisons to relatively coarse resolution simulations. • 

499 Our work here is ongoing [ Colarco et al., in preparation, 2006]. 

500 Qualitatively, the satellites (Figures 12 - 14) all show a similar aerosol feature over the 

501 same portion of southern Africa and eastern Atlantic Ocean (from approximately 0° to 30°S and 

502 12°W to 40°E). Over land, we find that the peak mean AOT observed by the satellites is further 

503 north and east (southwestern Zaire) relative to the base model (northeastern Angola). Over the 

504 ocean, we find that all the satellites detect higher monthly mean AOT values over the eastern 

505 Atlantic Ocean (approximately 5°S to 15°S and 0° to 15°E) relative to the base model. Despite 

506 differences in magnitude, the model agrees with the satellite retrievals in the general placement 

507 of the smoke plume. The hatched marks in Figures 12 - 14 represent locations with “No Data” 

508 (N.D.) or locations that did not pass the sampling requirements discussed above. MISR and EP- 

509 TOMS have more locations with no observations than MODIS which may be a result of the 

510 temporal and spatial coverage (i.e., larger pixel size) of each satellite. 

511 To quantitatively compare the satellite observations to the model. Figures 15 - 17 show 

512 scatterplots of the satellite versus model monthly mean AOT shown in Figures 12-14. Here we 

513 separate the comparison by land and ocean points. This division is particularly relevant for 

514 MODIS, which explicitly uses a different inversion algorithm over land than over ocean, but it 

515 also lets us compare the MISR and EP-TOMS performance over these two different surface 

516 classifications. The dashed lines represent the 1-to-l line and the solid lines represent the line of 

5 17 best fit. The correlation coefficient (R), the equation of the line of best fit, and the number of 

518 locations (listed in parentheses) are also reported in each scatterplot. For clarity, we grouped the 

519 uncertainty in the data into intervals of 0.2. Here, the error bars represent the average uncertainty 

520 in the data at each interval. In general, we find that our model results are linearly correlated to 
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521 the satellite observations, with stronger correlations over the ocean. Additionally, we find that 

522 the model has lower AOT values compared to all the satellite products, with larger 

523 disagreements in the ocean comparisons. Over land, we find that the model and MISR agree the 

524 best with the strongest correlations and slope nearest to 1 . However over the ocean, we find that 

525 the model agrees best to MODIS. We also find the largest disagreement between the model and 

526 EP-TOMS AOT over the ocean (i.e., weak correlation and significant offset in the slope of the 

527 line of best fit) which may be a result of the assumptions being made in the EP-TOMS retrievals. 

528 In particular, unfiltered sub-pixel clouds may be contributing to the larger AOT signal. 

529 There are at least three issues that could explain the model’s lower AOT values over the 

530 ocean. First, our sampling and aggregation strategy incorporated fewer pixels in a model grid 

531 box per day over the ocean than land due to sunglint (particularly in MODIS retrievals) and 

532 persistent cloud cover in certain regions, so we may be biased in that regard. Over land the 

533 average number of pixels in a model grid box per day during September 2000 was: 107 for 

534 MODIS, 48 for MISR, and 9 for EP-TOMS; while over the ocean: 84 for MODIS, 40 for MISR, 

535 8 for EP-TOMS. 

536 Second, the model may be deficient in that we are neglecting other important aerosol 

537 species, such as sea salt and dust. The mean Angstrom exponent at the 470/670 nm wavelength 

538 pair retrieved by MISR and MODIS during September 2000 over this region ranged between 0.2 

539 and 1.4. These values are commonly associated with large particles suggesting that dust and sea 

540 salt aerosols may be contributing to the higher AOT values, which would not be seen by our 

541 model since we are not simulating these aerosols. In addition, the model may not be properly 

542 accounting for significant microphysical evolution of smoke particles through, for example. 
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543 condensation of gases within the smoke plume, or the humidification of the smoke aerosols 

544 (discussed in more detail in section 3.2). 

545 A third point is that the high monthly mean AOT values detected by the satellite 

546 instruments over the Atlantic Ocean may be due to complications and/or errors in the satellite 

547 retrievals. A study completed by Ichoku et al. [2003] suggested that the global constant SSA of 

548 0.90 at visible wavelengths assumed in the smoke aerosol retrievals by MODIS is too high. 

549 Rather a SSA of 0.86, or smaller for fresh smoke aerosols, is more representative of biomass 

550 burning aerosols over southern Africa. Additionally, it is known that this region of the world has 

551 persistent, extensive stratus decks. The mean reflectivity from EP-TOMS for September 2000 

552 indicates high levels of reflectivity (~25 to 50%) over this region suggesting that unfiltered sub- 

553 pixel clouds may be contributing to the high AOT feature over the Atlantic Ocean. 

554 Since we observe inconsistencies in the comparisons between the model results and 

555 satellite observations, we compared MODIS to MISR (Figure 18). Figure 18 scatters the monthly 

556 mean AOT550 of MODIS versus MISR at the model resolution over the land and ocean. The 

557 monthly mean of the MISR and MODIS data was calculated by first aggregating the Level 2 

558 products to our model grid for each day using the same pixel thresholds discussed above. The 

559 monthly mean at each grid box was then constructed from the daily averaged grid box values 

560 where both MODIS and MISR made retrievals. The correlation coefficient (R), the equation of 

561 the line of best fit, and the number of grid points (listed in parentheses) are reported in each 

562 scatterplot. Again, the error bars represent the average uncertainty in the data at each interval of 

563 0.2. In general, we find that the datasets compare well over land and ocean with high 

564 correlations, low offsets and slopes near 1. However, MISR is lower than MODIS by 

565 approximately 25% over land and approximately 10% over the ocean. Abdou et al. [2005] found 
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566 similar comparisons between MISR and MODIS data collected during three months of 2002 

567 (March, June, and September). Since the satellite AOT observations are in good agreement over 

568 the ocean, the model is most likely missing information about the aerosols. However, the 

569 disagreement in the satellite AOT observations over land does not provide us with enough 

570 information to determine if the model is quantitatively doing a poor job. 

571 To better understand the model’s capability over land, we compared the simulated mean 

572 AOT to MISR, MODIS, and AERONET at the six study sites (Table 3). Note the satellite 

573 instruments should not be compared to AERONET quantitatively in this case since they have 

574 been degraded to our model resolution. Direct comparison of the satellite to AERONET is more 

575 correctly done using the higher resolution Level 2 products and examining the temporal 

576 correlation between the satellite and AERONET observations (e.g., as was done for MODIS in 

577 Ichoku et al. [2003]). Table 3 shows the simple mean of all days in the month where the MISR, 

578 MODIS, and AERONET datasets overlap, with AERONET and the model sampled as close to 

579 the satellite overpass time as possible. For AERONET, the mean represents the average of 

580 retrievals made within an hour of Terra’s overpasses time (i.e., between 800 UTC and 1000 

581 UTC), while the model average represents the results at 1200 UTC. Here we find that the model 

582 compares best to MISR at four of the six study sites (Etosha Pan, Mongu, Ndola, and Senanga) 

583 and compares best to MODIS at the remaining two sites. However, the model is consistently 

584 higher than MISR and lower than MODIS. Compared to AERONET, the model is lower at all 

585 sites except at Skukuza. Additionally, we find that MISR AOT is always lower than MODIS 

586 AOT. As noted above in section 3.1.1, we found that the optical properties at both Inhaca and 

587 Skukuza were different from the other sites since they are influenced by other aerosol sources 

588 (i.e., dust and sea salt). Therefore, we find it interesting that the model compares best to MODIS 
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589 at Inhaca and Skukuza, whereas the sites that are more clearly dominated by smoke compare 

590 better to MISR. Since the model behaves essentially the same at all sites, this discrepancy in the 

591 comparison of MOD IS and MISR probably reflects the differences in the retrieval techniques. 

592 This relationship may suggest that MISR is less sensitive to coarse mode particles, such as the 

593 aerosols found at Inhaca and Skukuza, while MODIS is not. These results may also suggest that 

594 the algorithms used by MISR may be more appropriate for retrievals that are dominated by a 

595 single aerosol source, where the MODIS algorithms may be more appropriate for retrievals with 

596 multiple aerosol sources. 

597 Since the EP-TOMS AOT is only provided at 380 nm, we constructed Table 4 to 

598 compare the model to the EP-TOMS and AERONET observations at the six study sites. Again to 

599 be consistent in the comparisons, the mean only includes days when all the instruments made 

600 retrievals, with AERONET and the model sampled as close to the satellite overpass time as 

601 possible. For AERONET, the mean represents the average of retrievals made within an hour of 

602 the Earth Probe’s overpass (i.e., between 930 UTC and 1 130 UTC), while the model mean 

603 represents the results at 1200 UTC. Here we find that the model mean AOT 380 compares best to 

604 AERONET at four of the six study sites (Etosha Pan, Inhaca, Senanga, and Skukuza). However 

605 the model has lower mean AOT 380 values at all sites except at Inhaca. We also find that the 

606 model is either higher or lower relative to the EP-TOMS mean values at the selected study sites. 

607 In this section we developed a method for quantitatively comparing the satellite retrieved 

608 AOT to the model. In general, we find that the model underestimates the AOT relative to all the 

609 satellite datasets. The model agrees best with MISR over land, MODIS over the ocean, while 

610 agrees the least with EP-TOMS AOT retrievals over both land and ocean. We also find that 

611 MODIS and MISR agree better over the ocean than over the land, despite them observing the 
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612 same scenes. These results may reflect the strengths and weaknesses of each satellite. A general 

613 problem found with all satellite sensors is sub-pixel cloud contamination, and is perhaps the 

614 worst with EP-TOMS because of its large pixel size. The better agreement between MODIS and 

615 MISR over the ocean suggests that the satellite algorithms use similar assumptions to retrieve the 

616 AOT. As Figure 18 indicates and noted in Table 3, MODIS and MISR do not agree with each 

617 over the land, and the model AOT compares better to MISR at Etosha Pan, Mongu, Ndola, and 

618 Senanga but better to MODIS at Skukuza and Inhaca. In the context of the AERONET 

619 observations we have already found that Skukuza and Inhaca are influenced by multiple aerosol 

620 sources and characterized by having more coarse mode aerosols than the other sites. This 

621 suggests that the various satellite retrieval algorithms have divergent capabilities in detecting the 

622 variability in the aerosol optical properties across the six land sites explored here. Over land sites 

623 dominated by fine mode aerosols, presumably smoke aerosols, we might reasonably expect more 

624 confidence in the MISR retrieval because of fewer assumptions on the behavior of the surface 

625 reflectance characteristics. However at land sites with multiple sources of coarse mode aerosols, 

626 most likely dust and sea salt aerosols, we might expect more confidence in MODIS because of 

627 its wider spectral range. In the end, the discrepancies between the satellites over land mean that it 

628 is unclear how well the model is really doing. Over the ocean, the better agreement between 

629 MODIS and MISR and the consistently smaller AOT in the model relative to those sensors 

630 suggests important contributions from other aerosol sources which are being neglected in the 

63 1 current version of the model. 

632 

633 3.1.4 Comparisons to EP-TOMS Aerosol Index 
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634 In this section, we compare the results from the base model to the aerosol index (AI) 

635 retrieved by EP-TOMS for September 2000 [ Torres et al. , 1998]. Briefly, the AI is a measure of 

636 the wavelength-dependent reduction of scattered radiance by aerosol absorption relative to a pure 

637 Rayleigh scattering atmosphere. The aerosol index is defined so that positive values generally 

638 correspond to UV-absorbing aerosols and negative values correspond to non-absorbing aerosols. 

639 However, when absorbing aerosols are at low altitudes (1.5 km or less), the contribution from 

640 absorption is masked by the relatively more dominant Rayleigh scattering, resulting in negative 

641 values of AI. Theoretically, the aerosol index has a nearly linear relationship to AOT (for 

642 AOT 3 so less than about 2) if the aerosols have a constant altitude and single scattering albedo 

643 [Torres et al., 1998, Hsu et al., 1999], Hsu et al. [1999] showed that the EP-TOMS AI 

644 measurements were linearly proportional to the AOT 380 derived independently from AERONET 

645 Sun-photometer measurements over regions of biomass burning at Mongu and Zambezi, Zambia, 

646 between July and October of 1996 and 1997. However, their relationship cannot be directly 

647 compared in our study because of changes in the definition of the TOMS Aerosol Index in the 

648 most recent version 8 of their algorithm [O. Torres, personal communication], resulting in higher 

649 AI values for a given AOT. For this reason, we determined a new relationship between AI and 

650 AOT 380 retrievals from AERONET made within one hour of the EP-TOMS overpass time at the 

651 study sites dominated by smoke in southern Africa (i.e., Etosha Pan, Mongu, Ndola, and 

652 Senanga) during September 2000. Here we found AI = 0.83xAOT38o- 0.4, with R = 0.66. 

653 In Figure 19 we scattered the monthly mean AI against the simulated monthly mean 

654 AOTsso- Here we restricted the geographic domain to a region dominated by biomass burning 

655 aerosols (i.e., 5°S to 20°S and 10°W to 35°E) and then separated the comparison by land (left 

656 column) and ocean (right column) points. The monthly mean AI was constructed using the Level 


30 



657 2 Aerosol Index product and same sampling method used to construct the monthly mean EP- 

658 TOMS AOT values. Here the model AOT 38 o was sampled near the satellite overpass time, and 

659 only where the satellite daily grid box averages passed the pixel threshold test. The correlation 

660 coefficient (R), the number of locations (listed in parentheses), and equation of the line of best fit 

661 are reported in each scatterplot. In general, we find that our model results are linearly correlated 

662 to the aerosol index, with stronger correlations over the ocean. The model’s regression line over 

663 land is not similar to that observed in the AERONET AOT 380 and EP-TOMS AI comparison 

664 discussed above (dashed line in Figure 19). The disagreement in the regression lines may be due 

665 to the model underestimating the AOT. As indicated above, when AERONET retrievals are 

666 greater than 1.0, the model underestimates the AOT by up to a factor of 2. The disagreement may 

667 also be due to the influence of other aerosol constituents that we are not simulating. Additionally, 

668 we find that EP-TOMS shows a larger AI over the ocean than over the land. Knowing that our 

669 simulations show smoke plumes rising in altitude as they move off the coast of Africa and that 

670 EP-TOMS observations are most sensitive to aerosols that are well above the surface, we would 

671 expect to see a larger signal by EP-TOMS off the coast, even though the column AOT may be . 

672 greatest over the continent. The strong signal over the ocean may also be a result of sub-pixel 

673 cloud contamination (discussed above) or the age of the smoke aerosols where the condensation 

674 of trace gases cause the aerosols to become more absorbing. 

675 

676 3.1.5 Comparisons to Airborne and Ground-Based Measurements of Particle Size Distributions 

677 In this section, we compare our simulated particle size distributions to AERONET 

678 retrievals and Passive Cavity Aerosol Spectrometer Probe 100X (PC ASP) measurements. The 

679 AERONET column-integrated aerosol volume size distributions are obtained from the Level 2 
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680 AERONET product which uses the application of the Dubovik and King [2000] algorithm to 

681 spectral almucantar sky radiances and spectral AOT from the CIMEL Sun-sky radiometers. The 

682 PC ASP instrument was onboard the C-130 of the UK Met Office airplane which flew within 

683 smoke source regions at altitudes between 208 m and 1000 m above ground level (AGL) during 

684 the SAFARI 2000 field campaign [Hay wood et al., 2003]. 

685 Figure 20 compares the daily mean simulated (solid line) and AERONET (dashed line) 

686 column-integrated volume size distributions at Etosha Pan on 13 September 2000. The model 

687 size distribution is normalized to the fine mode volume concentration of the AERONET 

688 retrieval. Compared to AERONET, we find that the model’s fine mode is broader and has a 

689 larger fine mode median radius. Additionally, the simulated volume distribution of particles with 

690 radius sizes greater than 2 (xm is up to a factor of 5 lower then AERONET. Knowing that Etosha 

691 Pan is influenced by sea salt and dust aerosols, which are large particles, the high coarse mode 

692 volume concentrations retrieved by AERONET may be a result of these aerosols. Since we are 

693 not simulating these aerosol sources, the model results will not be able to capture the high 

694 volume in the coarse mode. On the other hand, the base model is initialized with much higher 

695 volumes of coarse mode particles (Figure 2 - solid black line). This suggests that the model is 

696 either removing too many large particles or we are making incorrect assumptions about the 

697 smoke particle properties (i.e., density and shape), which is subsequently removing large 

698 particles too quickly. Further investigation of this issue is discussed in section 3.2. 

699 Figure 21 compares the daily mean volume size distribution simulated by the base model 

700 (solid line) to the volume size distribution fitted to parameters reported by Haywood et al. [2003] 

701 which are representative of average PC ASP measurements made in aged smoke plumes between 

702 208 m and 294 m AGL (dashed line) at Etosha Pan on 13 September 2000. The model size 
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703 distribution is normalized to the fine mode volume concentration of the PC ASP measurement. 

704 Again, we find that the model’s fine mode is broader and has more larger- sized fine mode 

705 particles than the PCASP measurement. However, the volume of particles with radius sizes 

706 greater than 1 pm is approximately a factor of 3 greater than PCASP, and the peak of the coarse 

707 mode simulated by the base model is shifted towards larger particles. Since the PCASP 

708 instrument was flown through relatively young smoke aerosols, the measurement may show 

709 larger volumes of particles at smaller radius sizes than the model. On the other hand, since the 

710 model was initialized with a particle size distribution representative of young and aged smoke, 

711 by this time the size distribution has had time to grow and coagulate resulting in larger volumes 

712 of particles at larger fine mode radius sizes relative to PCASP. Seeing that the model 

713 overestimates the coarse mode particle mass relative to the elevated altitude PCASP 

714 measurements but underestimates the coarse mode mass relative to the column-integrated 

715 AERONET retrieval suggests that the coarse mode particles observed exist primarily near the 

716 surface and may not be transported over long distances. 

717 

718 3.2 Results - Sensitivity Tests 

719 A series of sensitivity tests were conducted to understand the influence of various 

720 assumptions made in our model on the generality of our results. We partition these tests into two 

721 basic divisions. First, we consider a series of tests exploring the sensitivity of the model aerosol 

722 source to the timing and injection altitude of smoke emissions (Section 3.2.1). Second, we 

723 consider a series of tests exploring the sensitivity of the modeled aerosol properties to various 

724 assumptions about the particle microphysics, including initial particle size distribution, aging, 
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725 removal, and particle optics (Section 3.2.2). Each test is assigned a name, and a summary of the 

726 names and variances in each test are presented in Table 5. 

727 

728 3.2. 1 Sensitivity of Simulated Aerosol Distributions to Timing and Injection of Emissions 

729 In this section, we investigate how implementing a fluctuating daily emissions dataset, 

730 different aerosol emission injection altitudes, and diurnal cycle influences the model results. 

731 

732 3. 2.1.1 Sensitivity to Daily versus Monthly Aerosol Emissions 

733 An issue with using the GFED dataset is that it only provides monthly mean carbon 

734 emissions. An emission inventory with a finer temporal resolution may be important for 

735 realistically specifying the temporally varying simulated smoke optical properties. Therefore, 

736 rather than constantly emitting the same amount of aerosols throughout the month, we generated 

737 aerosol emissions that fluctuated day-by-day. The fluctuating daily aerosol emissions were 

738 constructed by correlating the emissions provided by the GFED to the Along Track Scanning 

739 Radiometer (ATSR) daily hot spot data. ATSR is a low-resolution sensor onboard the European 

740 Remote Sensing (ERS-2) satellite launched in April 1995 and has been used for fire location 

741 identification and environmental monitoring at visible and infrared wavelengths [Arino, 1995]. 

742 Here, we aggregated the daily ATSR hot spots to the same grid used in our model and emission 

743 source. The daily emissions for each grid cell are found by multiplying the monthly total 

744 emissions from the grid cell by the fraction of the monthly total fires occurring in the grid cell on 

745 that day. Potentially this underestimates the number of days on which our fires are active 

746 because heavy cloud or aerosol cover may mask fire hot spots that would otherwise be identified. 

747 On the other hand, the coarseness of the model grid might obviate this problem so that in an 
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748 average sense the day-to-day variability of fires occurring in the grid cell is about right. We note 

749 that the monthly mean emissions provided by the GFED dataset were constructed using the 

750 TRMM firecounts. However, we cannot use the TRMM firecounts here because of its irregular 

751 repeat cycle over a given region, which would cause it to sometimes miss a day entirely. In the 

752 monthly average sense, TRMM may be more useful because it sees more of the diurnal 

753 variability in emissions. We also note that we cannot use MODIS firecounts in this test because 

754 of problems with the MODIS detectors associated with hotspot detection during this time period 

755 [http://modland.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra]. 

756 In Figure 4 we compare the AOT 500 simulated in this daily emissions case (referred to as 

757 Daily Emissions and represented by the dashed line) to our base model and AERONET 

758 measurements. This comparison shows that fluctuating the aerosol emissions through the month 

759 does have some impact on the day-to-day variability of the aerosol load, particularly at Mongu, 

760 Ndola, and Skukuza. However in some cases, the AOT 50 o comparison in our Daily Emissions 

761 test is worse, suggesting that simply using hot spots to scale the gridded GFED emissions may 

762 miss important variables that determine the actual day-to-day emissions, such as the size of the 

763 fire, the vegetation burned, or the diurnal cycle of emissions. Furthermore, this case cannot be 

764 used to explain the missing aerosol in the model in mid-September at Mongu or Senanga, 

765 suggesting that there is either a missing source in the GFED or the model meteorology is 

766 somehow incorrect in this region at this time. Still, the Daily Emissions case and the base case 

767 are more similar than they are different, suggesting that the variation in the numbers of fires and 

768 timing of the fires throughout the month are not as important as the variability in the 

769 meteorological conditions for the local AOT 50 o- 

770 
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771 3. 2. 1.2 Sensitivity to the Diurnal Cycle of Aerosol Emissions 

772 As pointed out earlier, Eck et al. [2003] reported a diurnal variability in AOT retrieved by 

773 AERONET at several sites during SAFARI 2000. This may be a result of real diurnal patterns in 

774 the fires themselves, in which case the observations are seeing the plumes from nearby sources, 

775 or the diurnal variability may reflect the aerosol morphology or venting from the planetary 

776 boundary layer. In our base model where the particles were injected between 0600 UTC and 

777 1800 UTC every day, we found a very weak diurnal cycle in the simulated AOT at individual 

778 AERONET sites. Since we did not observe a strong diurnal trend in our base model, we 

779 conducted two additional tests where we varied the timing of emissions throughout the day. In 

780 one sensitivity test, we relaxed the base case diurnal cycle (referred to as “No DC” or no diurnal 

781 cycle) by evenly distributing the aerosol emissions throughout the day. In another test, we 

\ 

782 injected all the aerosols at the 1500 UTC timestep on each day (referred to as “1500 UTC”). 

783 Relative to the base case, changing the diurnal cycle of the aerosol emissions only had a small 

784 effect on the model extinction profile at the lowest altitudes (between 0 and 0.5 km) and almost 

785 no effect above that (not shown). In the test where nocturnal emissions were allowed, this 

786 suggests that those emissions can stay near the surface at night, but are quickly mixed throughout 

787 the boundary layer during the day. Overall, we are unable to reproduce the diurnal variability in 

788 AOT observed by Eck et al. [2003]. It is not clear that changes in aerosol morphology through 

789 any diurnal variations in the optics are to blame, as our computed aerosol optical properties agree 

790 reasonably well with the observations available (see Section 3.2.5). Therefore, it seems more 

791 likely that at the coarse spatial resolution of our model and the smoke emissions from point 

792 sources that AERONET is sensitive to are dispersed quickly in our grid. 

793 
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794 3. 2. 1.3 Sensitivity to Aerosol Injection Height 

795 To test the model’s sensitivity to the altitude at which the aerosols are initially distributed 

796 vertically above their source regions, we ran two additional simulations that were based on 

797 information provided from the literature. During SAFARI 2000, trace gases and aerosols over 

798 the subcontinent were consistently observed as trapped stable layers at approximately 850 hPa 

799 (~1.5 km), 700 hPa (-3.0 km), and 500 hPa (-5.5 km) within the lower troposphere [Stein et al . , 

800 2003, Schmid et al., 2003, McGill et al., 2003]. The injection altitude of aerosol emissions was 

801 assumed to be associated with buoyancy generated by the release of heat in the fires, and the 

802 amount of heat release was found to be dependent on the amount and type of vegetation being 

803 burned, and whether the combustion was smoldering or flaming [ Andreae and Merlet, 2001]. 

804 Using this information, we ran one simulation with aerosol emissions centered at three model 

805 layers and referred to as the 3-Layers model (Layers included: 1 .42 km (0.41 km thick), 3.06 km 

’ x 

806 (0.67 km thick), and 5.64 km (0.98 km thick)). In another simulation, we emitted the aerosols 

807 within one model layer, centered at 2.42 km (0.58 km thick) and this test is referred to as the 2 

808 km-Layer model. 

809 Figures 10 and 1 1 show mean vertical extinction profiles for the test simulations. For a 

810 one-day comparison (Figure 10), we find that the overall vertical profiles of the aerosol are 

811 almost u n i f orm with altitude. The results from the 2 km-Layer model is not shown in Figure 10 

812 since the results are nearly identical to the base model. However in the 3-Layers model, the 

813 aerosols get carried downwind more rapidly due to the higher wind speeds aloft, thereby 

814 producing extinction values that are too low relative to the observations and base model. In 

815 Figure 1 1, we compare the tests to the monthly mean MPLNET lidar observations at Skukuza. 

816 Here a different picture emerges where we find that the simulated vertical profiles from our two 
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817 elevated injection tests are not similar to the base model. Instead, the aerosol appears to reside at 

818 an elevated altitude relative to both the base case and the observations. Our base model of 

819 uniform mixing throughout the depth of the PBL is therefore the most comparable to the 

820 observations. 

821 

822 3.2.2 Sensitivity of Simulated Aerosol Properties to Microphysical Processes 

823 In this section we discuss the implications and limitations of various assumptions made 

824 about the aerosol microphysics in our model, including choices of initial particle size 

825 distribution, resolution of bin sizes, optical properties, and aerosol aging and removal processes. 

826 

827 3.2.2. 1 Sensitivity to Treatment of Particle Size Distribution 

828 CARMA is a bin-resolving cloud and aerosol microphysics module. Incorporation of 

829 such a module into a large-scale transport-modeling framework can easily become 

830 computationally prohibitive. To reduce this computational expense, CARMA can be adjusted to 

83 1 use the fewest number of size bins and microphysical processes necessary to carry out a given 

832 problem. In contrast, a more typical chemical transport model undertaking the problem we 

833 presented here might have a greatly reduced aerosol microphysical modeling capability 

834 compared to CARMA. As an example, the NASA Global Ozone Chemistry Aerosol Radiation 

835 and Transport (GOCART) model [Chin et al. 2002] neglects the explicit treatment of aerosol 

836 particle size distributions for carbonaceous particles, instead partitions the aerosols into 

837 hydrophobic and hydrophilic organic carbon (OC) and black carbon (BC) components. The 

838 subdivision of the aerosol mass to hydrophilic and hydrophobic components affects the aerosol 

839 optics, as well as the aerosol lifetime since only the hydrophilic fraction is scavenged by 
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840 precipitation. The time-scale for conversion of hydrophobic to hydrophilic aerosols in these 

841 models is typically short (approximately 1 day) and is intended to represent condensation of 

842 gases and water absorption onto the particles, but does not allow for explicit treatment of other 

843 aging mechanisms such as coagulation or the treatment of the aerosol optics as internal mixtures. 

844 On the other hand, it should be pointed out that GOCART is more complete than our model in 

845 the sense that it treats all major tropospheric aerosols species. 

846 The principle strength of our model is the explicit treatment of specific microphysical 

847 aging processes through coagulation and sedimentation, which is a step towards evolving 

848 internally mixed aerosol species. A question we present here is how many size bins are required 

849 to adequately constrain our problem. Based on earlier simulations (e.g., Colarco et al. [2004]), 

850 we chose to use 16 size bins spaced from 0.01 - 10 Jim radius in our base model. As a limiting 

851 case, we also considered a case with 22 size bins spaced between 0.05 - 15 pm radius which is 

852 identical to the bin separation AERONET uses to report their size distributions. Using the 

853 refractive index and particle size distribution retrieved on 16 September 2000 at Ndola, we 

854 compute the Angstrom exponent for the 440 and 870 nm wavelength pair ( 0 ( 440 / 870 ) and curvature 

855 term at 380, 500, and 870 nm wavelengths (a’ 380 / 500 / 870 ) for the 16-bin and 22-bin cases, and then 

856 compared the values to the AERONET observation on this particular day. Here, we find that the 

857 CX 440/870 is 1-93, 1.97, and 1.96, respectively, for the AERONET observation, 16-bin case, and 22- 

858 bin case. Additionally, the curvature terms are 1.21, 1.03, and 1.12, respectively, for the same 

859 three cases. Here, the computed Angstrom exponents are similar to the observed value, and 

860 presumably the differences can mainly be attributed to actual particle non-sphericity reflected in 

861 the measured values. On the other hand, compared to the observations, the computed curvature 

862 terms are approximately 20% and 10% smaller in the 16-bin case and 22-bin case, respectively. 
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863 Note that the curvature term was defined in terms of the 380, 500, and 870 nm wavelengths, and 

864 of those three wavelengths AERONET only retrieves refractive indices at 870 nm. Here, we 

865 interpolate the 500 nm refractive index from the 440 and 670 nm retrieved values. Additionally, 

866 the 380 nm refractive index is assumed to be the same as the 440 nm retrieved value, which 

867 seems to be a reasonable assumption based on the discussion in Torres et al. [2005]. The results 

868 of our calculations are essentially the same when we employ the wavelength independent 

869 refractive index used in our previous base model calculations. 

870 The discrepancy in the computed curvature terms above can be attributed to the number 

871 of size bins used to resolve the fine mode portion of the particle size distribution. We note that 

872 there is essentially no contribution to either the Angstrom exponent or curvature term from 

873 particles smaller than 0.05 jim radius or larger than 10 jam radius. The fine mode particle 

\ 

874 scattering area is primarily between particles 0.05 and 0.6 |im radius, which comprises six size 

875 bins in our 16-bin case and nine size bins in our 22-bin case. From sensitivity tests with more 

876 bins (not shown), we find that increasing the number of size bins in the fine mode improves the 

877 accuracy of the computed optics. Here we find a convergence in the accuracy of the computed 

878 optical properties when ten or more bins are used in the fine mode particle size distribution (i.e., 

879 between 0.05 and 0.7 jam). The improved accuracy of the optics by adding more size bins is also 

880 illustrated in the comparison of the initial particle size distribution simulated Angstrom exponent 

881 for the 380 and 440 nm wavelength pair, which is 1.48, 1.52, and 1.47 for the observation, 16-bin 

882 case, and 22-bin case, respectively. The discussion above focuses on the uncertainty in the 

883 prescribed initial particle size distribution used in the model, and the model’s ability to reproduce 

884 the observed optical properties. In the following discussion, we consider different assumptions 

885 about the initial particle size distribution and how that evolves in the model. 
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886 , Our simplest sensitivity test is called the Mass model. Here we simply integrate the mass 

887 across all the size bins in our base model and compute the aerosol optical properties following 

888 Chin et al. [2002], assuming the aerosols were completely hydrophilic because of the short 

889 timescale for conversion. We also consider different assumptions about the initial particle size 

890 distribution based on retrievals by AERONET (column-integrated) or measurements made by the 

891 airborne PC ASP instrument {in situ) [Haywood et al., 2003] over smoke source regions during 

892 the SAFARI 2000 field campaign. In the PCASP model we used a particle size distribution that 

893 was a fit to parameters determined from mean PCASP observations taken within a smoke plume 

894 at Otavi, Na mi bia on 13 September 2000. The AERONET-Mean model used a particle size 

895 distribution that was fitted to the mean parameters of smoke dominated days retrieved at Ndola, 

896 Zambia by AERONET in September 2000. The 22-Bins model resolves the particle size 

897 distribution using the AERONET reported bins as described above. Table 5 lists the model 

898 names, input parameters that were changed in each model, and the number fractions, geometric 

899 standard deviations and median number radii associated with each particle size distribution. 

900 Figure 2 presents the column-integrated volume size distributions that were implemented into the 

901 models. 

902 Table 2 lists the monthly mean AOT 500 , CC 380 / 440 , a ’> SSA 440 , SSA 670 , SSA 870 , and SSA 1020 

903 of smoke dominated days for selected sensitivity tests at Ndola and Inhaca. Note that not all of 

904 the results are shown for Inhaca since they show the same conclusions as at Ndola. In general, 

905 we find that the simulated mean AOT500 and SSA are the most sensitive to Mass and PCASP 

906 parameterizations. Compared to the base model, the Mass model mean AOT 500 is approximately 

907 50% higher at both study sites. The reason for this is that the Mass model implicitly assumes that 

908 all the aerosol mass is in the fine mode, which has greater aerosol extinction per unit mass than 
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909 the coarse mode. However, the low monthly mean Angstrom exponent (a = 1.27) of the Mass • 

910 model at both locations suggests that its assumed fine mode particle size distribution has a 

911 slightly larger median radius than the base model or observations. We also find that the Mass 

912 model is approximately 55% more absorbing than the base model at both study sites. 

913 We also find that the PCASP model monthly mean AOT 500 is approximately 28% lower 

914 than the base model at Ndola. The high mean Angstrom exponent (a = 1.66) of the PCASP 

915 model suggests that this model is transporting more smaller sized particles relative to the base 

916 model, but the lower mean curvature value suggests a greater relative contribution of coarse 

917 mode particles to the total AOT. Figure 22 compares the simulated column-integrated volume 

918 size distributions of the base and PCASP models to the AERONET retrieval on 16 September 

919 2000 at Ndola. The model size distributions are normalized to the fine mode volume 

\ 

920 concentration retrieved by AERONET. Here we find that the PCASP model has a smaller fine 

921 mode median radius compared to base model, but a larger volume (mass) concentration at radius 

922 sizes between 1.0 and 10 pm, which explains the higher Angstrom exponent (smaller fine mode) 

923 but lower optical depth and curvature term (more of the mass in the coarse mode). We also find 

924 that the PCASP model is roughly 25% more absorbing then the base model at both study sites. 

925 The monthly mean AOT500, curvature term and SSA results from the AERONET-Mean 

926 model were not significantly different from the base model results at Ndola (see Table 2). 

927 However, the mean Angstrom exponent [380/440] of the AERONET-Mean model is higher than 

928 the base model suggesting that the AERONET-Mean model has a smaller fine mode median 

929 radius. 

930 As discussed above, we find that the 22-Bins model does impact our computed optical 

931 properties. Here we find that the monthly mean a 38 o /440 is slightly lower and the curvature value 
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932 is higher .than the base model which has 16 size bins. These results are consistent with our 

933 findings of the initial calculations of the 16-bin case and 22-Bin case discussed above. We also 

934 find that the 22-Bins model monthly mean CC 380/440 and curvature values differ from the 

935 AERONET observations. By increasing the bin resolution of the model, we are able to observe 

936 more details about the simulated particle size distribution than illustrated in the base model 

937 alone. Here, the 0 : 330/440 of the 22-Bins model and AERONET observation is 1.32 and 1.43, 

938 respectively, which suggests that the particle fine mode median radius is too large in the model 

939 relative to the observations. The curvature term of the 22-Bins model and AERONET 

940 observation is 1.38 and 1.23, respectively, suggesting that relative to AERONET too much of the 

941 aerosol optical thickness is due to the fine mode. This point is qualitatively bom out in almost all 

942 of our simulations because the model does not maintain the observed concentration of coarse 

943 mode particles (see Figure 22 for the 16-bin base case and discussion in Section 3. 2.2. 3). To 

944 explain the lower Angstrom exponent we recall that the initial particle size distribution already 

945 included some effects of aging and coagulation of smoke particles so that further coagulation has 

946 made the fine mode median radius slightly too large. 

947 Overall, these results suggest that the transport of a size distribution and the appropriate 

948 bin resolution of the accumulation mode in our model are important in order to obtain optical 

949 properties that are comparable to AERONET. The Mass model produced mean AOT 50 o values 

950 that were higher and aerosols that were more absorbing than both the base model and 

95 1 AERONET at the study sites. Additionally, the Mass model has the disadvantage that it cannot 

952 capture the aging of the particles as they evolve over several days. Westphal et al. [1991] and 

953 Colarco et al. [2004] simulated the evolution of smoke clouds over several days and found that 

954 the aerosol optical properties changed as the particles evolved by coagulation. Westphal et al. 
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955 [1991] demonstrated that small changes in the geometric mean radius yielded large changes in - 

956 extinction and absorption efficiencies and therefore in the specific extinction and absorption. 

957 However, the sign of the change depended on the size distribution. Eck et al. [1999] found that 

958 the mode size was positively correlated with the aerosol optical thickness. Reid et al. [ 1998] 

959 suggested that this increase in radius was due to the aging processes of smoke aerosols, whereby 

960 the particle size increases with time because of both coagulation and gas-to-particle conversion. 

961 Our model does not include gas to particle conversion, but the impact of coagulation is discussed 

962 below. 

963 

964 5.2. 2.2 Sensitivity to Microphysical Processes 

965 The size of the aerosol is important when calculating optical properties because different 

966 particle sizes scatter different amounts of radiation. Several observational studies have shown 

967 that aerosol size distributions shift from fine mode particle sizes to accumulation mode particle 

968 sizes as the aerosol ages due to both coagulation and gas-to-particle conversion [ Reid and Hobbs, 

969 1998; Reid et al., 1998; Eck et al, 2001; Eck et al., 2003; Haywood et al., 2003]. Simulations of 

970 smoke aerosols over multi-day periods have also suggested that coagulation is necessary in order 

971 to explain the change in particle size distributions and optical properties as the smoke ages 

972 [ Westphal et al., 1991; Colarco et al., 2004]. In this section, we investigate the model’s 

973 sensitivity to coagulation by conducting a simulation where the particles were not allowed to 

974 coagulate (referred to as the No Coagulation model). 

975 In our No Coagulation simulation, we find that the mean AOT 500 of smoke dominated 

976 days is approximately 10% lower than the base model at Ndola. Additionally, the Angstrom 

977 exponent is higher while the curvature term is lower relative to both the base model and the 
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978 AERONET observations (Table 2). The computed SSA for the No Coagulation model was 

979 slightly lower (i.e., more absorbing) than both the base model and the AERONET observation, 

980 although both the base and No Coagulation models are within the uncertainty of the AERONET 

981 retrieval error bars. The interpretation of these results is rather complicated by the uncertainty in 

982 the resolution of the spectral optical properties in our base model. At first glance, the higher 

983 Angstrom exponent and lower curvature values of the No Coagulation model suggest too small 

984 of a fine mode median radius and excessive contribution from coarse mode particles to the total 

985 AOT. However, initializing our model with a particle size distribution that already includes some 

986 aged aerosols and then allowing the particles to coagulate further exaggerates the effects of 

987 coagulation making it difficult to reproduce the correct optical properties. Therefore, the No 

988 Coagulation model where we already implicitly include coagulation may be a more reasonable 

989 approximation for the aerosol optical properties over southern Africa. Here the base model is 

990 probably more reasonable for long-range transport of the aerosols over longer periods of time, 

991 which is consistent with Westphal et al. [1991] and Colarco et al. [2004], 

992 

993 3. 2.2.3 Sensitivity to Removal Schemes 

994 The properties of smoke particles vary depending on fuel type and moisture, combustion 

995 phase, and wind conditions. These properties can also change rapidly as they disperse. In 

996 general, smoke particles are composed of -60% organic carbon and -5-10% black carbon, 80- 

997 90% of their volume is in the accumulation mode, reported density varies between 0.79 and 1.35 

'i 

998 g cm , and studies have shown smoke particles to have a variety of morphologies such as chain 

999 aggregates, solid irregulars, and more liquid/spherical shapes [ Reid et al., 2004]. About -10% of 
1000 the smoke volume consists of a variety of coarse mode particles (-2.5-15 jim). These coarse 
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1001 mode ash particles (2 < d p < 20 (im) have been observed more than fifty kilometers from the - 

1002 source [ Reid et al., 2004], However, the density and exact shape of these giant ash particles is 

1003 unknown to our knowledge. 

1004 Since the model’s volume concentrations of particles greater than 2 jim were found to be 

1005 approximately a factor of 5 too low, this suggests that we are most likely not treating the large 

1006 particles correctly in our model. There are two possibilities that may explain this model error. 

1007 First, given the large uncertainty in the coarse mode of the AERONET size distribution 

1008 retrievals, we may not be representing the coarse mode particles properly in our model. Second, 

1009 our removal processes (i.e., dry deposition routine) may not be appropriate for this region. 

1010 Currently, we are treating all of the smoke particles by applying the same particle density (1.35 g 

101 1 cm -3 ) and spherical shape to all the particle sizes. However the coarse mode particles could 

\ 

1012 consist of ash particles, which could have low density and irregular shape. For these reasons, we 

1013 conducted two additional simulations; one in which we changed the particle density to 0.675 g 

1014 cm" 3 (i.e., half the base model density and referred to as Density model) and another in which we 

1015 made the particles flat-plates by following Fuchs [1964] (referred to as Shape model). On the 

1016 other hand, we may be using the appropriate smoke properties in our base model, however the 

1017 dry deposition parameterization used to calculate the dry deposition velocities may be removing 

1018 the large particles in our model too vigorously. Therefore, we conducted two additional 

1019 simulations; one where we did not allow the particles to be removed via dry deposition (referred 

1020 to as No Dry Dep. model) and another where we implemented an alternative method to calculate 

1021 the dry deposition velocities (referred to as Zhang et al. [2001] model). The alternative dry 

1022 deposition calculation we examined in this study follows the parameterization described by 
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Zhang et.pl. [2001] which is more robust than the current dry deposition routine by Shao [2000]. 
In particular, the Zhang et al. [2001] considers surface type where Shao [2000] does not. 

Figure 23 shows the normalized daily mean column-integrated volume size distribution 
retrieved by AERONET and simulated by the test models at Ndola, Zambia on 16 September 
2000. The model distributions are normalized to the fine volume concentration retrieved by 
AERONET. For clarity, we did not present the results from the Shape or Zhang et al. [2001] 
models since the fine mode results were similar to the base model and the coarse mode results 
were similar to the Density model. In general, we find that different particle shapes and densities, 
and an alternative method in computing the dry deposition velocities improve the volume of the 
coarse mode particles but not enough to be comparable to the AERONET retrievals. The 
combination of all three modifications (density, shape, and dry deposition method) in a single 
simulation provided results that were the most comparable to AERONET (not shown but similar 
to the Density fine mode and the No Dry Dep. coarse mode). Of the results shown in Figure 23, 
the only method that reproduced similar volume concentrations of particles greater than 2 pm 
relative to AERONET was the No Dry Dep. model. Since we are able to maintain the large 
particles in the No Dry Dep. model, this suggests that our model has problems transporting large 
particles. On the other hand, since the initial particle size distribution in the model may represent 
both local and transported aerosols, the coarse mode of the transported aerosol may have been 
even more pronounced at the source. Thus, our base case may not be initialized with enough 
coarse mode particles. Altogether, without more information about the properties of the coarse 
mode smoke particles, in particular about the ash particles, we do not know which method is the 
most appropriate to use in our model. Further investigation of the smoke size distribution is 
subject of future research. 
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Of the tests conducted in this section, we found that the simulated optical properties were 
significantly impacted by the change in particle density. Since we want to maintain the same 
injected total particle mass, decreasing the particle density by half doubles the number of 
particles injected which results in increasing both the rate of coagulation of fine mode particles 
and the computed AOT, but decreasing the computed Angstrom exponent. The optical properties 
of the other tests in this section are essentially the same as the base model. The slightly lower 
curvature terms of these tests relative to the base model are expected and due to the additional 
coarse mode particles each test provides. 

3.2.2A Sensitivity to Smoke Optics 

In this section, we investigate how the choice of refractive index and relative humidity 
influences the model results. 

To examine the model’s sensitivity to choice of refractive index, we tested four different 
refractive indices that are representative of biomass burning aerosols (Table 5). The refractive 
indices tested here were either calculated by Haywood et al. [2003] or retrieved by AERONET. 
Haywood et al. [2003] calculated a refractive index of 1.54 - 0.025/ at 550 nm for young smoke 
aerosols. This refractive index was applied to all wavelengths (wavelength independent (WI)) 
and referred to as the WI-Haywood et al. [2003] model. We also selected a refractive index 
retrieved on 16 September 2000 at Ndola by AERONET (N re f = 1.52 - 0.019/ at 500 nm). Again, 
we chose this particular day since it had the highest AOT 500 which suggests that this day was 
largely dominated by smoke aerosols. This refractive index was applied to all wavelengths and 
referred to as the WI-AERONET model. We also explored wavelength dependent (WD) 
refractive indices, which are listed in Table 5. The WD-AERONET model applied refractive 
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1069 indices retrieved by AERONET on 16 September, and the WD-AERONET-Mean model applied 

1070 the mean refractive indices of smoke dominated days retrieved by AERONET at Ndola, Zambia. 

1071 Here we find that the monthly mean AOT 500 and Angstrom exponent [380/440] of smoke 

1072 dominated days, are not significantly sensitive to the choice of complex refractive index (Table 

1073 2). However, we did find that the simulated mean curvature values and SSA were sensitive to our 

1074 choice of refractive index (see Table 2 and Figure 24). Compared to the base model, the WI 

1075 models have higher curvature values, while the WD models have lower curvature values. 

1076 Additionally, all but the WI - Haywood et al. [2003] model have higher mean SSA (i.e., more 

1077 absorbing aerosols) than the base model. Overall, we find that the WI refractive indices are the 

1078 most comparable to the AERONET observations where the WD cases are not. 

1079 To test the model’s sensitivity to relative humidity, we considered a test where we 

1080 assumed that the relative humidity at all locations was zero (referred to as No RH model). 

1081 Compared to the base model, we find that when we do not treat the humidification of aerosols, 

1082 the mean AOT 50 o is lower and the particles become more absorbing at Ndola. However, the 

1083 Angstrom exponent and curvature terms remain relatively unchanged (Table 2). 

1084 Overall, we find that WI refractive indices and treating relative humidity produces optical 

1085 properties that are the most comparable to AERONET observations at Ndola. 

1086 

1087 4. Conclusions 

1088 In this study, we investigated the transport and optical properties of biomass burning 

1089 aerosols using a three-dimensional aerosol microphysical, transport, and radiation model. Our 

1090 first goal in this study was to determine if using estimated aerosol emissions, measured particle 

1091 size distributions, and estimated injection altitudes collected during the SAFARI 2000 field 
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campaign as input parameters in our model allowed us to reproduce the observed optical 
properties. Our second goal was to better understand the microphysics needed to accurately 
simulate biomass burning aerosol properties over our domain. The principle findings of our study 
are summarized as follows: 

1. Comparisons to AERONET AOT: We find that our base model, which uses a 
satellite-derived monthly inventory of biomass burning emissions and AERONET 
derived aerosol optical properties, was sufficient to reproduce the majority of the day- 
to-day variability in AOT at the six AERONET sites used in this study. In certain 
instances where the model significantly underestimated the aerosol loading with 
respect to AERONET, we found from satellite imagery that the model had misplaced 

the smoke plume relative to observations, suggesting errors in the model dynamics. 

\ 

Sensitivity tests in which we varied the diurnal emissions of the smoke or the timing 
of emissions throughout the month using daily satellite firecounts showed little effect 
on the daily variability of simulated AOT at the selected AERONET locations. 
Although these results do not exclude possible errors in the biomass burning emission 
inventory used, they suggest that the aerosol loading over the continent is most 
strongly controlled by dynamics, which is consistent with results from Myhre et al. 
[2003]. On a monthly basis, the model reproduces mean AOT values which are lower 
than those observed at the six AERONET sites, suggesting either an uncertainty in the 
magnitude of our aerosol emissions or that we are neglecting contributions from other 
aerosol species. 
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2. Comparisons to Satellite Observations ofAOT : Comparisons of the model AOT to 
satellite observations showed that our model is capable of simulating smoke plume 
features similar to those observed by SeaWiFS, MODIS, MISR, and EP-TOMS over 
southern Africa and the Atlantic Ocean. Qualitatively, this suggests that the input 
parameters and microphysical processes in our model are adequate for this study 
region. In order to quantitatively compare our model to the satellite observations we 
constructed monthly mean regional AOT maps in which the model observations were 
weighted and sampled in manner consistent with the available satellite observations. 
We find that the simulated monthly mean AOT is linearly correlated to the satellite 
observations over the land and ocean, but that the model is consistently lower than the 
satellite observations. Over land, we find that our model results compare best to 
MISR with the strongest correlations and slope nearest to 1 . At the AERONET sites, 
the model compares best to MISR at the four sites dominated by biomass burning, 
while the model compares best to MODIS at the two remaining sites which are 
influenced by multiple aerosol sources. Over the ocean, the model compares the best 
to MODIS. Over both land and ocean we find the largest disagreements between the 
model and EP-TOMS. These results may reflect the differences in the retrieval 
technique of each satellite. Here the MISR algorithms may be more appropriate for 
retrievals dominated by a single aerosol source, while the MODIS algorithms may be 
more appropriate for retrievals with multiple aerosol sources. Due to the large pixel 
size of the EP-TOMS retrievals, the large disagreement between the model and EP- 
TOMS may be a result of sub-pixel clouds that are not being removed from their 
retrievals. While sub-pixel cloud contamination may affect all the aerosol retrievals. 
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this might be expected to the most pronounced over the ocean where we find 
generally fewer pixels available for constructing our monthly means than over land. 
The disagreement among the model and satellite observations may also be due to the 
contributions of other aerosol species to the column AOT, which might be more 
pronounced over the ocean than land. Additionally, we cannot exclude the errors in 
the assumptions of the aerosol optical properties used in the satellite retrievals or in 
the model. In particular, MODIS uses a different retrieval algorithm over land points 
than over the ocean. Overall, the consistently lower AOT simulated by the model 
relative to MODIS and MISR over the ocean, and the good agreement between 
MODIS and MISR over the ocean suggests that the model is missing important 

aerosols or processes. However, the disagreement between MODIS and MISR over 

\. 

the land does not provide us with enough information to determine the model’s 
performance over continental Africa. 

/ 

3. Comparisons to Airborne and Ground-Based Measurements of Extinction: The 
simulated vertical extinction profiles from our base model compare favorably to the 
AATS-14, MPL, and CPL measurements. The minor differences found between the 
simulations and observations on a single day may be a result of the model’s coarse 
resolution and source function. Given that our source function cannot represent the 
likely variations in daily smoke emissions, it is challenging to reproduce individual 
days. However, the model’s vertical resolution is sufficient in representing long time- 
scale averaged vertical profiles. When exploring simulations with elevated aerosol 


52 



1159 

1160 
1161 
1162 

1163 

1164 

1165 

1166 

1167 

1168 

1169 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 
1181 


injections, we found that our best results were with our base model assumption of 
injections uniformly mixed within the planetary boundary layer. 

4. Comparisons to Measurements and Retrievals of Particle Microphysical and Optical 
Properties: While the base model compares favorably to daily and monthly 
observations of AOT by AERONET, it does not capture the daily variability in the 
Angstrom exponent and retrieved single scatter albedo. It is likely that the model 
remains consistent because of its relatively invariant aerosol composition. At the 
selected sites, the observed fluctuations may be due to the diumally varying 
contributions from dust and other aerosol species that are not simulated in this study. 
On a monthly mean basis, the modeled spectral dependence of the AOT and single 
scatter albedo compare favorably with the AERONET observations at Ndola (and the 
other smoke do m inated sites) and slightly less at Inhaca (which has important 
contributions from non-smoke aerosols). Comparisons of the simulated particle size 
distribution to AERONET retrievals and airborne PCASP measurements reveal that 
the model base case evolves a slightly larger fine mode median radius than both sets 
of observations. Additionally, these particle size distribution comparisons show that 
the model overestimates the coarse mode particle mass relative to the elevated 
altitude PCASP measurements but underestimates the coarse mode mass relative to 
the AERONET retrieval. Assuming that there is no error in the AERONET volume 
distribution of the coarse mode, the latter point suggests that the coarse mode 
particles observed exist primarily near the surface and may not be transported over 
long distances. We also found that the model is evidently too aggressive in removing 
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these coarse mode particles, and sensitivity tests suggest that different assumptions . 
about the composition, shape, or density of the coarse mode particles would help with 
this issue. On the other hand, the larger fine mode median radius in the base model 
may be due to our assumption about the initial particle size distribution, which 
already includes some effects of the smoke aging. These general points are further 
illustrated in the optical calculations presented in Table 2. The results of these tests 
indicate that the simulated optical properties are relatively insensitive to choices of 
particle refractive index or humidification (i.e., composition), but are sensitive to the 
choice of initial particle size distribution and, as discussed in Section 3.2.2. 1, to the 
number of size bins used to resolve the fine mode fraction of the particle size 
distribution. Here we find that initializing the model with the airborne PCASP 
measured particle size distribution does not represent the evolution of the coarse 
mode, and the Mass model does not reproduce the observed optical properties. 
Initializing the model with the AERONET size distribution retrieved at Ndola leads to 
too large of a fine mode median radius relative to the observations. Overall, for a 
model that treats coagulation, the implication is that the fine mode needs to be 
initialized to a smaller median radius than what is reported by AERONET. Sensitivity 
tests of the optical calculations further suggest that at least ten size bins are needed to 
adequately resolve the optical properties of the evolving fine mode. 


Since Africa experiences the most extensive biomass burning in the world, this work is a 
step toward improved quantification of smoke optical properties for use in calculating the direct 
effects of these aerosols on climate, and suggests methods for incorporating the particle size 
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1205 distribution of smoke aerosols into microphysical models attempting to simulate the aging of 

1206 smoke particles. However, we did find inconsistencies in the AOT land comparisons among 

1207 AERONET, MISR, MODIS, EP-TOMS data that need to be better understood in order for us to 

1208 feel confident in using them to aid us in modeling and reducing the uncertainty in aerosol 

1209 forcing. These inconsistencies and further investigation of aging mechanisms inspire further 

1210 research. 
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1454 September 2000. 

1455 Table 4. Presents the number of days included in the mean AOT and the mean AOT observed by 

1456 AERONET, EP-TOMS, and simulated by the base model at the six study sites for September 

1457 2000. 

1458 Table 5. Sensitivity test names and modifications made relative to Base Model. 
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1459 Figure Captions 

1460 Figure la. Satellite image of the “River of Smoke” on 4 September 2000. This image shows an 

1461 example of biomass burning aerosols being transported across southern Africa and over the 

1462 Indian Ocean [Swap et al. [2003], SeaWiFS Project, NASA/Goddard Space Flight Center, and 

1463 ORB IMAGE], 

1464 Figure lb. Simulated daily mean AOT at 550 nm On 4 September 2000. 

1465 Figure 2. Vertically column-integrated volume size distributions used in the simulations 

1466 conducted in this study. The solid black line represents the base model particle size distribution 

1467 fitted to parameters from an AERONET almucantar scan at Ndola, Zambia on 16 September 

1468 2000. The gray dashed line represents the PC ASP model size distribution fitted to parameters 

1469 from mean PCASP measurements at Otavi, Namibia on 13 September 2000. The gray solid line 

1470 represents the AERONET -Mean model size distribution fitted to parameters from the mean 

1471 AERONET almucantar scans at Ndola, Zambia on smoke dominated days in September 2000. 

1472 The gray dashed line represents the 22-Bins model size distribution. The parameters for each of 

1473 the models are listed in Tables 1 and 4. 

1474 Figure 3. Mean interpolated aerosol emissions using GFED dataset for September 2000 and 

1475 locations where the data used in our paper was collected. 

1476 Figure 4. Daily mean AOT 50 o simulated by the base model (solid line) and Daily Emissions 

1477 model (dashed line), and all AOT 500 retrievals from AERONET (+) at each study site for 

1478 September 2000. 

1479 Figure 5. Scatterplot of daily mean AOT 500 retrieved by AERONET versus simulated daily 

1480 mean AOT 500 for smoke dominated days at the six study sites for September 2000. The dashed 

1481 line represents the 1-to-l line and the solid line represents the line of best fit through all points. 
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1482 The equation of the line of best fit and correlation coefficient (R) is reported in the figure. Each 

1483 study site is represented by a different symbol: Etosha Pan (square), Inhaca (triangle), Mongu 

1484 (diamond), Ndola (plus sign), Senanga (asterisk), and Skukuza (cross). 

1485 Figure 6 . Scatterplot of daily mean SSA 440 retrieved by AERONET versus simulated daily mean 

1486 SSA 440 for smoke dominated days at the six study sites for September 2000. The dashed line 

1487 represents the 1-to-l line and the solid line represents the line of best fit through all points. The 

1488 equation of the line of best fit and correlation coefficient (R) is reported in the figure. Each study 

1489 site is represented by a different symbol: Etosha Pan (square), Inhaca (triangle), Mongu 

1490 (diamond), Ndola (plus sign), Senanga (asterisk), and Skukuza (cross). The single bar at the top- 

1491 left represents the error in the AERONET retrievals of +0.03. 

1492 Figure 7. Scatterplot of daily mean Angstrom exponent [440/870 nm] retrieved by AERONET 

1493 versus simulated daily mean Angstrom exponent [440/870 nm] for smoke dominated days at the 

1494 six study sites for September 2000. The dashed line represents the 1-to-l line and the solid line 

1495 represents the line of best fit through all points. The equation of the line of best fit and 

1496 correlation coefficient (R) is reported in the figure. Each study site is represented by a. different 

1497 symbol: Etosha Pan (square), Inhaca (triangle), Mongu (diamond), Ndola (plus sign), Senanga 

1498 (asterisk), and Skukuza (cross). 

1499 Figure 8 . Monthly mean spectral dependence of AOT of smoke dominated days for September 

1500 2000 at Inhaca and Ndola. The solid lines represent the results from the base model and the 

1501 asterisks connected with a dashed line represent the AERONET retrievals. 

1502 Figure 9. Monthly mean spectral dependence of SSA of smoke dominated days for September 

1503 2000 at Inhaca and Ndola. The solid lines represent the results from the base model and the 

1504 asterisks connected with a dashed line represent the AERONET retrievals. 
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1505 Figure 10. Mean extinction profiles at Mongu and Senanga on 6 September 2000. The simulated 

1506 profiles represent the average of the results at 0600 and 1200 timesteps, the CPL profiles 

1507 represents the average of four samples between 0800 and 1100 UTC (black line), and the AATS 

1508 profiles represent one sample around 0900 UTC (red line). The blue line represents the 

1509 simulated results from the base model, and the green line represents the simulated results from 

1510 the 3 -Layers model. 

1511 Figure 11. Mean extinction profiles at Skukuza for September 2000. Mean includes: 1-3, 6 , 7, 9, 

1512 10, 13, 14, 17 September 2000. Base model is represented by the solid black line, 3-Layers 

1513 model is represented by the solid red line, 2 km-Layer model is represented by the solid blue 

1514 line, and MPL retrievals are represented by the dashed line. 

1515 Figure 12. Comparisons of monthly mean AOT 55 o observed by MODIS (left) and simulated by 

v 

1516 our base model (right) for September 2000 over southern Africa. The hatch marks represent 

1517 locations with no data (N.D.). 

1518 Figure 13. Comparisons of monthly mean AOT 550 observed by MISR (left) and simulated by our 

1519 base model (right) for September 2000 over southern Africa. The hatch marks represent locations 

1520 with no data (N.D.). 

1521 Figure 14. Comparisons of monthly mean AOT 3 so observed by EP-TOMS (left) and simulated by 

1522 our base model (right) for September 2000 over southern Africa. The hatch marks represent 

1523 locations with no data (N.D.). 

1524 Figure 15. Scatterplots of monthly mean AOT550 observed by MODIS versus model. The domain 

1525 is divided into to two regions: land (left) and ocean (right). The dashed line represents the 1-to-l 

1526 line and the solid line represents the line of best fit. 
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1527 Figure 16. Scatterplots of monthly mean AOT 550 observed by MISR versus model. The domain 

1528 is divided into to two regions: land (left) and ocean (right). The dashed line represents the 1-to-l 

1529 line and the solid line represents the line of best fit. 

1530 Figure 17. Scatterplots of monthly mean AGT 380 observed by EP-TOMS versus model. The 

1531 domain is divided into to two regions: land (left) and ocean (right). The dashed line represents 

1532 the 1-to-l line and the solid line represents the line of best fit. 

1533 Figure 18. Scatterplots of monthly mean AOT 55 o observed by MODIS versus MISR at points 

1534 over land (left) and ocean (right) for September 2000. The dashed line represents the 1-to-l line 

1535 and the solid line represents the line of best fit. 

1536 Figure 19. Scatterplots of monthly mean EP-TOMS AI versus simulated monthly mean AOT 380 

1537 for September 2000. The geographic domain is restricted to a region dominated by biomass 

1538 burning aerosols (5°S to 20°S and 10°W to 35°E) and then separated into land (left) and ocean 

1539 (right) points. The solid lines in both figures represent the line of best fit. The dashed line in land 

1540 plot represents the regression line from the comparison of AERONET AOT 38 o and EP-TOMS AI 

1541 over smoke dominated locations in September 2000. 

1542 Figure 20. Daily mean column-integrated volume size distribution simulated by the base model 

1543 (solid line) and the daily mean column-integrated volume size distribution fitted to parameters 

1544 from AERONET almucantar scans (dashed line) at Etosha Pan on 13 September 2000. The 

1545 model size distribution is normalized to the fine volume concentration of the AERONET 

1546 retrieval. 

1547 Figure 21. Daily mean volume size distribution simulated by the base model (solid line) and the 

1548 volume size distribution fitted to parameters reported by Haywood et al. [2003] which are 

1549 representative of average PCASP measurements within an aged smoke (dashed line) at Etosha 
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1550 Pan on 13 September 2000. The model size distribution is normalized to the fine volume 

1551 concentration of the PC ASP measurement. 

1552 Figure 22. Normalized column-integrated daily mean volume size distributions at Ndola, Zambia 

1553 on 16 September 2000. The model size distributions are normalized to the fine volume 

1554 concentration of the AERONET retrieval. The black solid line represents the base model, gray 

1555 solid line represents the PCASP model, and the dashed black line represents the AERONET 

1556 retrieval. 

1557 Figure 23. Normalized daily mean column-integrated AERONET and simulated volume size 

1558 distributions at Ndola, Zambia on 16 September 2000. The model size distributions are 

1559 normalized to the fine volume concentration of the AERONET retrieval. The solid black line 

1560 represents the base model, gray line represents Density model, the gray dashed line represents 

\ 

1561 the No Dry Dep. model, and the dashed black line represents the AERONET retrieval. 

1562 Figure 24. AERONET and simulated mean spectral dependence of SSA of smoke dominated 

1563 days for September 2000 at Inhaca and Ndola. The solid line represents the results from the Base 

1564 model, the dotted line represents the WD- AERONET model, the dashed-dotted line represents 

1565 the WI-AERONET model, the dashed line represents the WD-AERONET-Mean model, and the 

1566 asterisks connected with dashed line represent the AERONET retrievals. 
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1568 


Table 1. Input Parameters used in Base Model. 


Microphysical Processes 

Advection 
Coagulation 
Dry Deposition 
Wet Removal 
Sedimentation 

Emission Dataset 

GFED [van der Werfet al., 2003] 

Emission Factor: 10 g TPM/kg of DM 

Emission Rate over Entire 
Model Domain 

Constant throughout month 

August: 2.05e-7 g cm -2 s - ' 
September: 3.64e~7 g cm -2 s’" 1 

Injection Height 

Mixed between Surface & PBLH 

\ 

Peak height of the MATCH PBL ranged between 1 
km and 3.7 km at the study six sites 

Diurnal Cycle 

Aerosol emitted for 12-hour period during daylight hours 

Particle Density 

1.35 gem -3 

Initial Size Distribution 

Radius Range [pm] 

# Bins 

f 

Mode 

N f 

r N [pm] 

c 

0.01 to 10 

16 

1 

0.99986 

0.079 

14 

2 

0.00014 

0.695 

2.] 

Dependence of Size on 
Relative Humidity 

Follows parameterization described by Magi & Hobbs [2003] 

Refractive Index 

Wavelength independent 

N ref = 1.51-0.024/ 
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Table 2. Monthly mean AOT 500 , ocpsomo], a’ pso/soara], SSA440, SSA 670 , SSA 87 o, and SSA 102 o of smoke dominatec 
days from various simulations and AERONET at Inhaca and Ndola for September 2000. 









































































































































1571 Table 3. Presents the number of days included in the mean AOT and the mean AOT observed by 

1572 AERONET, MISR, MODIS, and simulated by the base model at the six study sites for September 200 0. 


Site 

# of Days 

Model r550] 

MISR [550] 

MODIS [550] 

AERONET [500] 

Etosha Pan 

7 

0.38 

0.37 

0.68 

0.51 

Inhaca 

6 

0.36 

0.28 

0.37 

0.38 

Mongu 

7 

0.59 

0.54 

0.65 

0.79 

Ndola 

6 

0.49 

0.47 

0.62 

0.80 

Senanga 

7 

0.60 

0.54 

0.73 

0.75 

Skukuza 

5 

0.46 

0.37 

0.47 

0.30 


1573 

1574 Table 4. Presents the number of days included in the mean AOT and the mean AOT observed by 

1575 AERONET, EP-TOMS, and simulated by the base model at the six study sites for September 2000. 


Site 

# of Days 

Model [380] 

EP-TOMS [380] 

AERONET [380] 

Etosha Pan 

23 

0.50 

0.71 

0.54 

Inhaca 

11 

0.61 

0.51 

0.58 

Mongu 

20 

1.05 

1.18 

1.35 

Ndola 

18 

1.02 

1.11 

1.15 

Senanga 

18 

1.07 

0.83 

1.16 

Skukuza 

11 

0.67 

0.53 

0.69 
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Table 5. Sensitivity test names and modifications made relative to Base Model. 


Model Name 

Parameter Changed in Base Model 

New Parameter 

PC ASP 

Size Distribution (Mode/Number 
Fraction/Radius[N]/Std. Dev.) 

1 

0.9957 

0.10 pm 

1.; 

2 

0.0041 

0.22 pm 

l.f 

3 

0.0002 

1 .00 pm 

1.9 

AERONET - Mean 

Size Distribution (Mode/Number 
Fraction/Radius[N]/Std. Dev.) 

1 

0.9999 

0.08 pm 

1 .i 

2 

0.0001 

0.73 pm 

2.1 

Mass 

Size Distribution 

Transport mass following Chin et al. [2001] 

22 - Bins 

Radius Range 

0.05 to 15 pm 

Number of Size Bins 

22 Size Bins 

3 Layers 

Injection Height 

Aerosols emitted in three model layers (1. 5/3/5 km 

2 km-Layer 

Injection Height 

Aerosols emitted in one model layer (2 km) 

No DC 

Diurnal Cycle 

Aerosols emitted evenly throughout the day 

1500 UTC 

Diurnal Cycle 

Aerosols emitted at 1 500 UTC time step each day 

Daily Emissions 

Emission Rate 

GFED emissions linked to ATSR Fire Count Data 

No Coagulation 

Microphysical Process 

No Coagulation 

No RH 

Relative Humidity 

Relative Humidity assumed to be zero 

WI - Haywood et al. T20031 

Refractive Index (A/N„ f ) 

Wavelength Independent 

1 .54-0.0251 

WI- AERONET 

Refractive Index (X/N ref ) 

Wavelength Independent 

1.52-0.019i 

. 

WD - AERONET- Mean 

' 

Refractive Index (A/N ref ) 

380 nm 

1.51-0.0241 

440 nm 

1.51-0.0241 

670 nm 

1.52-0.0201 

870 nm 

1.53-0.0201 

1020 nm 

1.54-0.01 61 

WD- AERONET 

Refractive Index (?jN ref ) 

380 nm 

1.52-0.0191 

440 nm 

1.52-0.0191 

670 nm 

1.52-0.0161 

870 nm 

1.53-0.0171 

1020 nm 

1.53-0.0161 

Density 

Density of Particle 

Density: 0.675 g/cm 3 

Shape 

Shape of Particle 

Particle Shape: Flat Plates 

No Dry Dep. 

Microphysical Process 

No dry deposition in model 

Zhang et al. [2001] 

Microphysical Process 

Dry deposition routine from Zhang et al. [2001] 



1579 Figure la. Satellite image of the “River of Smoke” on 4 September 2000. This image shows an 

1580 example of biomass burning aerosols being transported across southern Africa and over the 

1581 Indian Ocean [Swap et al. [2003], SeaWiFS Project, NASA/Goddard Space Flight Center, and 

1582 ORBIMAGE]. 
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1600 Figure lb. Simulated daily mean AOT at 550 nm on 4 September 2000. 
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1602 
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1612 
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1615 


Figure 2. Vertically column-integrated volume size distributions used in the simulations 
conducted in this study. The solid black line represents the base model particle size distribution 
fitted to parameters from an AERONET almucantar scan at Ndola, Zambia on 16 September 
2000. The gray dashed line represents the PCASP model size distribution fitted to parameters 
from mean PCASP measurements at Otavi, Namibia on 13 September 2000. The gray solid line 
represents the AERONET-Mean model size distribution fitted to parameters from the mean 
AERONET almucantar scans at Ndola, Zambia on smoke dominated days in September 2000. 
The gray dashed line represents the 22-Bins model size distribution. The parameters for each of 
the models are listed in Tables 1 and 4. 
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Figure 3. Mean interpolated aerosol emissions using GFED dataset for September 2000 and 
locations where the data used in our paper was collected. 
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1622 


Figure 4. Daily mean AOT 500 simulated by the base model (solid line) and Daily Emissions 
model (dashed line), and all AOT 500 retrievals from AERONET (+) at each study site for 
September 2000. 


1.8 

1.2 
0.6 
0 

1.8 

1.2 
0.6 
0 

1.8 

1.2 
0.6 
0 

5 10 15 20 25 30 5 10 15 20 25 30 

Day in September 2000 




79 





1623 

1624 

1625 

1626 

1627 

1628 


1629 

1630 

1631 

1632 

1633 

1634 

1635 

1636 

1637 

1638 


1639 

1640 


Figure 5. Scatterplot of daily mean AOT 50 o retrieved by AERONET versus simulated daily 
mean AOT 500 for smoke dominated days at the six study sites for September 2000. The dashed 
line represents the 1-to-l line and the solid line represents the line of best fit through all points. 
The equation of the line of best fit and correlation coefficient (R) is reported in the figure. Each 
study site is represented by a different symbol: Etosha Pan (square), Inhaca (triangle), Mongu 
(diamond), Ndola (plus sign), Senanga (asterisk), and Skukuza (cross). 



Figure 6 . Scatterplot of daily mean SSA 440 retrieved by AERONET versus simulated daily mean 
SSA440 for smoke dominated days at the six study sites for September 2000. The dashed line 
represents the 1-to-l line and the solid line represents the line of best fit through all points. The 
equation of the line of best fit and correlation coefficient (R) is reported in the figure. Each study 
site is represented by a different symbol: Etosha Pan (square), Inhaca (triangle), Mongu 
(diamond), Ndola (plus sign), Senanga (asterisk), and Skukuza (cross). The single bar at die top- 
left represents the error in the AERONET retrievals of +0.03. 
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1641 Figure 7. Scatterplot of daily mean Angstrom exponent [440/870 nm] retrieved by AERONET 

1642 versus simulated daily mean Angstrom exponent [440/870 nm] for smoke dominated days at the 

1643 six study sites for September 2000. The dashed line represents the 1 -to- 1 line and the solid line 

1644 represents the line of best fit through all points. The equation of the line of best fit and 

1645 correlation coefficient (R) is reported in the figure. Each study site is represented by a different 

1646 symbol: Etosha Pan (square), Inhaca (triangle), Mongu (diamond), Ndola (plus sign), Senanga 

1647 (asterisk), and Skukuza (cross). 
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Figure 8. Monthly mean spectral dependence of AOT of smoke dominated days for September 
2000 at Inhaca and Ndola. The solid lines represent the results from the base model and the 
asterisks connected with a dashed line represent the AERONET retrievals. 
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1655 Figure 9. Monthly mean spectral dependence of SSA of smoke dominated days for September , 

1656 2000 at Inhaca and Ndola. The solid lines represent the results from the base model and the 

1657 asterisks connected with a dashed line represent the AERONET retrievals. 
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1661 Figure 10. Mean extinction profiles at Mongu and Senanga on 6 September 2000. The simulated 

1662 profiles represent the average of the results at 0600 and 1200 timesteps, the CPL profiles 

1663 represents the average of four samples between 0800 and 1100 UTC (black line), and the AATS 

1664 profiles represent one sample around 0900 UTC (red line). The blue line represents the 

1665 simulated results from the base model, and the green line represents the simulated results from 

1666 the 3-Layers model. 
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Figure 11. Mean extinction profiles at Skukuza for September 2000. Mean includes: 1-3, 6, 7, 9, 
10, 13, 14, 17 September 2000. Base model is represented by the solid black line, 3-Layers 
model is represented by the solid red line, 2 km-Layer model is represented by the solid blue 
line, and MPL retrievals are represented by the dashed line. 



Extinction [525 nm] (km' 1 ) 


83 



1675 

1676 

1677 


1678 

1679 

1680 
1681 
1682 


1683 

1684 

1685 

1686 
1687 


1688 


Figure 12. Comparisons of monthly mean AOT 550 observed by MOD IS (left) and simulated by 
our base model (right) for September 2000 over southern Africa. The hatch marks represent 
locations with no data (N.D.). 
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Figure 13. Comparisons of monthly mean AOT 550 observed by MISR (left) and simulated by our 
base model (right) for September 2000 over southern Africa. The hatch marks represent locations 
with no data (N.D.). 
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Figure 14. Comparisons of monthly mean AOT 380 observed by EP-TOMS (left) and simulated by 
our base model (right) for September 2000 over southern Africa. The hatch marks represent 
locations with no data (N.D.). 
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1689 Figure 15. Scatterplots of monthly mean AOT 550 observed by MODIS versus model. The domain 

1690 is divided into to two regions: land (left) and ocean (right). The dashed line represents the 1-to-l 

1691 line and the solid line represents the line of best fit. 
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Figure 16. Scatterplots of monthly mean AOT 550 observed by MISR versus model. The domain 
is divided into to two regions: land (left) and ocean (right). The dashed line represents the 1-to-l 
line and the solid line represents the line of best fit. 
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1699 Figure 17. Scatterplots of monthly mean AOT 380 observed by EP-TOMS versus model. The 

1700 domain is divided into to two regions: land (left) and ocean (right). The dashed, line represents 

1701 the 1-to-l line and the solid line represents the line of best fit. 
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1703 Figure 18. Scatterplots of monthly mean AOT 550 observed by MODIS versus MISR at points 

1704 over land (left) and ocean (right) for September 2000. The dashed line represents the 1 -to- 1 line 

1705 and the solid line represents the line of best fit. 
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Figure 19. Scatterplots of monthly mean EP-TOMS AI versus simulated monthly mean AOT 380 
for September 2000. The geographic domain is restricted to a region dominated by biomass 
burning aerosols (5°S to 20°S and 10°W to 35°E) and then separated into land (left) and ocean 
(right) points. The solid lines in both figures represent the line of best fit. The dashed line in land 
plot represents the regression line from the comparison of AERONET AOT 380 and EP-TOMS AI 
over smoke dominated locations in September 2000. 
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Figure 20. Daily mean column-integrated volume size distribution simulated by the base model 
(solid line) and the daily mean column-integrated volume size distribution fitted to parameters 
from AERONET almucantar scans (dashed line) at Etosha Pan on 13 September 2000. The 
model size distribution is normalized to the fine volume concentration of the AERONET 
retrieval. 



Figure 21. Daily mean volume size distribution simulated by the base model (solid line) and the 
volume size distribution fitted to parameters reported by Haywood et al. [2003] which are 
representative of average PCASP measurements within an aged smoke (dashed line) at Etosha 
Pan on 13 September 2000. The model size distribution is normalized to the fine volume 
concentration of the PCASP measurement. 
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Figure 22. Normalized column-integrated daily mean volume size distributions at Ndola, Zambia 
on 16 September 2000. The model size distributions are normalized to the fine volume , 
concentration of the AERONET retrieval. The black solid line represents the base model, gray 
solid line represents the PCASP model, and the dashed black line represents the AERONET 
retrieval. 



Figure 23. Normalized daily mean column-integrated AERONET and simulated volume size 
distributions at Ndola, Zambia on 16 September 2000. The model size distributions are 
normalized to the fine volume concentration of the AERONET retrieval. The solid black line 
represents the base model, gray line represents Density model, the gray dashed line represents 
the No Dry Dep. model, and the dashed black line represents the AERONET retrieval. 
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Figure 24. AERONET and simulated mean spectral dependence of SSA of smoke dominated 
days for September 2000 at Inhaca and Ndola. The solid line represents the results from the Base 
model, the dotted line represents the WD-AERONET model, the dashed-dotted line represents 
the WI- AERONET model, the dashed line represents the WD-AERONET-Mean model, and the 
asterisks connected with dashed line represent the AERONET retrievals. 
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